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ABSTRACT 

This  work  investigates  the  effects  of  digital  processing 
in  radiometers.   It  deals  mainly  with  two  digital  versions 
of  a  total  power  radiometer.   The  first  consists  of  RF, 
Mixer  and  IP  sections  followed  by  an  analog  to  digital 
converter.   All  further  processing  is  done  in  a  digital 
computer.   The  second  version  consists  of  RF,  Mixer  and  IF 
sections  followed  by  a  square  lav;  detector,  RC  filter  and 
analog  to  digital  converter.   From  this  point  on  the  pro- 
cessing is  done  by  a  digital  computer.   A  figure  of  merit 
is  defined  based  on  the  performance  of  an  analog  total  power 
radiometer.   Exact  results  are  obtained  for  the  figure  of 
merit  of  the  first  digital  version.   For  the  second,  an 
approximate  solution  is  obtained.   The  effects  of  saturation 
and  finite  step  size  of  the  quantizer  were  taken  into  con- 
sideration for  the  above  results.   The  performance  of  digital 
balanced-Dicke  and  noise-adding  radiometers  is  investigated 
using  the  above  results.   The  effects  of  digital  filtering 
on  the  performance  of  a  radiometer  is  considered. 
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I.   INTRODUCTION 

A  radiometer  is  basically  a  highly  sensitive  and  stable 
noise  receiver.   The  principle  of  operation  lies  in  the  fact 
that  any  object  at  temperatures  above  absolute  zero  radiates 
energy  in  the  form  of  electromagnetic  waves.   Until  some 
years  ago,  radiometers  were  used  mainly  as  radio  telescope 
receivers  and  constituted  the  main  tool  of  Radio  Astronomy 
[Ref.  1].   Lately,  their  use  has  expanded  to  many  other 
areas.   The  measurement  of  atmospheric  temperature  [Refs.  2 
and  3] »  detection  of  air  turbulence  in  the  Troposphere 
[Ref.  4],  airborne  mapping  [Ref.  5],  the  measurement  of 
absolute  radiation  from  a  projectile  flow  field  [Ref.  6], 
and  passive  detection  [Ref.  7]  can  be  cited  as  examples. 

The  type  of  processing  involved  in  a  radiometer  is 
especially .fitted  to  digital  methods.   A  number  of  insti- 
tutions have  integrated  the  radiometer  with  the  digital 
computer  in  their  applications.   The  purpose  of  this  work 
is  to  investigate  the  effects  of  this  integration,  and  to 
provide  guidelines  to  the  design  of  digital  radiometers. 

The  main  body  of  the  investigation  will  deal  with  two 
digital  versions  of  a  total  power  radiometer  (TPR)  shown  in 
Figure  1.   The  first  one  (Figure  2)  consists  of  RF,  Mixer, 
and  IF  sections,  followed  by  an  Analog  to  Digital  Converter 
(ADC)  section.   All  further  processing  is  done  in  a  digital 
computer.   The  second  one  (Figure  3)  consists  of  RF,  Mixer, 


and  IF  sections  followed  by  a  square  law  detector,  a  low 
pass  filter,  and  an  ADC.   From  that  point  on,  the  processing 
is  done  by  a  digital  computer. 

The  main  difference  of  the  two  cases  is  the  location  of 
the  ADC.   The  second  version  (Figure  3)  is  presently  being 
used  and  its  main  advantages  over  the  first  version  (Figure 
2)  is  that  it  requires  a  much  lower  sampling  rate. 

The  results  obtained  for  a  total  power  digital  radiometer 
can  be  used  to  determine  performance  of  other  types  of 
radiometers.   Performance  of  digital  Dicke  and  noise-adding 
radiometers  will  be  considered. 
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II.   PERFORMANCE  FIGURE  OF  A  RADIOMETER. 
DEGRADATION  FACTOR 

By  definition,  the  sensitivity,  or  minimum  detectable 
temperature,  of  a  radiometer  is  taken  as  that  change  in 
operating  temperature  T  that  causes  a  change  in  the  expected 
value  of  the  output  equal  to  one  standard  deviation  of  the 
output  [Ref.  1].   Refer  to  Figure  4.   It  shows  the  expected 
value  and  the  standard  deviation  of  the  output  voltage  as 
functions  of  temperature  for  a  typical  radiometer. 


outpo  t" 
vol+ctoe. 


ECU 


FIGURE  4 
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Prom  Figure  4,  It  follows  that  (to  first  order) 


aeci]  .  ffim 


AT 


(1) 


op 


where  E[«]  =  Expected  value  of  [•],  and  T   is  operating 
temperature.   But  by  definition 


AE[I]  =  a- 


=  dE[I] 

T       dT 
op 


op 


(2) 


Therefore 


AT  = 


"I 

dELIJ 
dT 


'op 


(3) 


A  performance  figure  can  be  defined  as, 


AT 

T 
op 


T 

dE[I] 
dT 


op 


(4) 


It  can  be  shown  [Ref.  1]  that  for  an  analog  total  power 
radiometer  (Figure  1),  using  the  fact  that  the  inverse  of 
the  integration  time  is  much  smaller  than  the  bandwidth  of 
the  RF,  Mixer,  IF  filter,  that 


E[I]  =  TkT     /  |H  (f)|   df 


°P  _co"   '  ° 


(5) 


aT2   =  2xk2T  2   /  |HQ(f)|4  df 


(6) 
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where 

t  =  Integration  time 

H  (f)  =  Band-pass  filter  (RF,Mix,IF)  frequency 
response 

k  =  Boltzmann1 s  constant 
Substituting  into  (4) 

\  AT   _    1 


Top    (Bet)^2 


(7) 


where 


(   /  |Ho(f)|2  df)2 

be  =  — °= ; <8> 

/  |H  (f)|?  df 
0 

The  performance  figure  of  the  analog  TPR,  (7),  will  be 
used  as  a  basis  for  evaluating  the  performance  of  a  digital 
TPR.   Hence,  the  Degradation  Factor  can  be  defined  as 


-,  _  Performance  figure  of  Digital  Radiometer 
Performance  figure  of  Analog  Radiometer 


Intuitively,  F  should  be  equal  or  greater  than  one  since 
as  the  sampling  rate  of  the  ADC  goes  to  infinity  and  the 
step  size  of  the  quantizer  goes  to  zero,  both  digital 
radiometers  are  reduced  to  an  analog  TPR. 

Some  comments  are  appropriate  at  this  point.   First, 
(5)  and  (6)  show  that  for  analog  radiometers  E[I]  and  a_ 


1A 


are  linear  functions  of  temperature  and  go  through  the 
origin.   Hence  equation  (4)  reduces  to 


iT  -  ^T  (9) 


T 
op 


The  addition  of  a  quantizer,  due  to  saturation,  makes 
E[I]  and  a_  nonlinear  with  respect  to  temperature,  and  the 
performance  figure  is  thus  temperature  dependent.   There- 
fore (4)  must  be  used  instead  of  the  better-known  (7). 
Secondly,  the  parameter  T   in  equation  (4)  is  not  a 
convenient  choice,  since  for  a  given  T   ,  different  gain 
settings  of  the  RP,  Mixer,  IP  filter  will  change  the 
performance  of  the  radiometer.   Referring  to  Figures  2  and 
3  it  is  easily  shown  that 


T  =  a  aY  2  (Figure  2)        (10) 

c 


T=  3  az  (Figure  3)        (11) 


where  a  and  3  are  proportionality  constants.   Using  the 

2 
parameters  a    or  a  as  appropriate,  instead  of  T,  norma- 

c 

lizes  the  problem  of  gain  setting  out  of  the  performance 

figure . 
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Then,  for  the  radiometer  of  Figure  2 


AT 
^op 


dE[IJ 


da. 


aX  2  at  T 
c      op 


(12) 


and  for  the  radiometer  of  Figure  3 


AT 

r 

op 


dELl] 
da„ 


a7  at  T   , 
Z     op 


(13) 
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III.   DEGRADATION  FACTOR  CALCULATION  WITH  NO  QUANTIZER 

An  ADC  performs  two  basic  operations  on  a  signal.   One 

is  sampling,  and  the  other  is  amplitude  quantization.   In 

this  chapter  the  effects  of  sampling  alone  are  considered. 

Hence  it  is  assumed  that  each  sample  can  be  represented 

without  any  error. 

\ 
A.   IF  SAMPLING 

Refer  to  Figure  2.   Assume  that  at  the  antenna  white 

gaussian  noise  is  present  with  double-sided  power  spectral 

density  No/2  watts/hertz.   Then  X(t)  is  a  Narrow  Band 

Gaussian  Process  with  known  power  spectral  density,  and 

can  be  represented  as 


X(t)  =  X  (t)cos27rf  t  -  X  (t)sin2TTf  t         (14) 

C  OS  o 


where  X  (t)  and  X  (t)  are  baseband  Gaussian  processes  and 
s         c 

f  is  the  IF  center  frequency  of  the  radiometer.   It  can 

be  shown  [Ref.  8]  that  if  X(t),  X  (t)  and  X  (t)  are  to  be 

c         s 

wide  sense  stationary,  then 


RX  (y)  =  RX  (y)  (15) 

c        s 


RY(y)    =   RY    (u)cos2Trf  y   +   RY  v    (y)sin2TTf  u  (16) 

c  s   c 
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where  R  (y)  is  the  autocorrelation  function  of  the  process 

X(t)  and  RY  y.    (y)  is  the  crosscorrelation  function  of  the 

s  c 

'processes  X  (t)  and  X  (t). 
c         s 

Furthermore  if  H  (f)  is  symmetrical  about  f  (this  is 

o  o 

the  only  case  considered  in  this  work) ,  then 


Rx  x  (y)  =  0  (17) 

c  s 


and  RY  (y)  is  given  by 
Xc 


RY  (y)  =  2  /  S  (f)cos[2Tr(f-f  )y]df     (18) 
Ac       0   x  ° 


and  Sv  (f),  the  power  spectral  density  of  X  (t),  is  the 

A  C 

c 

Fourier  Transform  of  equation  (18). 

Physically,  for  the  case  of  interest  (i.e.  X(t),  X  (t) 

and  X  (t)  wide  sense  stationary,  and  H  (f)  symmetrical  about 
s  o 

f  ),  S„  (f)  is  Just  twice  the  low  pass  version  of  SY(f). 

O      A  A 

C 

Referring  to  Figure  2,  the  low  pass  filters  filter  out  the 
double-frequency  terms  coming  out  of  the  multipliers.   Then 


I  =  I±   +  I2  (19) 


I-,  -  }■     2  X„2(kA)  (20) 

x   n  k=1  c 


since  X  (kA)  =  q  (kA)  and 
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*o  =  h        l      X  2(kA)  (21) 

d       n  k=l  s 

since  X  (kA)  =  q  (kA). 
s        s 


Using  the  fact  that  for  Gaussian  random  variables 

E[X1X2X3X4]  =  E[X1X2]EOC3Xi|]+E[X1X3]E[X2X4]+E[X1Xi|]E[X2X3] 

(22) 
and  expanding  the  double  summation,  it  follows  that 

2oY   ^     n-1 
aT  d   =  — 2_  +  1     e   (1-£)RY  d(kA)         (23) 
Xl      n     n  k=l    n  Xc 


and 


EC^]  =  ax  2  (24) 


For  the  random  variable  I?  the  same  results  are  obtained 

since  Rx  (u)  =  Rx  (y) 
c        s 


Then 


E[II]=  2aY  2  (25) 

2     Xc 


and 


aj2  =  oI  2  +  aI  2  +  2E[I1I2]  -  2E[I1]E[I2] 

1      2  "  (26) 
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Expanding  E[I  I?]  and  using  the  fact  that  R„  „  (y)  =  0, 

c  s 
the  last  two  terms  of  (26)  cancel  out. 

Using  equations  (25)  and  (26)  in  (12)  it  follows  that 


op         naY  k=l  c 

c 

or 


1   (^-)2  =  |  +  -^  ^   (1  -  |)Rx2(kA)       (28) 
op         xax  k=l  c 

c 

where 


t  =  nA  (29) 


The  degradation  factor  squared  is 


5         2B„A  n-1  .  .    9 

F^  =  B  A  +  — \     I      (1  -  £A)R  ^(kA)         (30) 

ax   k=l  c 
c 

RY  (y)  is  a  function  of  H  (f),  therefore  (30)  will  be  used 

Ac  ° 

with  four  different  frequency  characteristics. 


1. ■  Ideal  Bandpass  Filter 
Refer  to  Figure  5 


N  B  sirnrBy 


V(p>  =  -^m —  (31) 
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and 


V  =  NQ  B  (32) 


Substituting  (3D  and  (32)  into  (30) 


F2  =  BpA  +  2BT,AnE1(l-^)  ^llLI^  (33) 

E      E  k=l    T   (irkAB)2 


But 


BE  =  B  =  W  (34) 

where  BE  is  the  equivalent  bandwidth  from  (8)  and  W  is  the 


half- power  bandwidth. 

Then  (33)  can  be  rewritten  as 


t?2   TTA  ,  otta  v   sin  irkW    2(WA)    r   ksin  TrkW    r->c\ 
F  =  WA  +  2WA  E   5-  -      '        l      — p—   (35) 

k=l   (TrkWA)^     TW   k=l   (TTkWA) 


Since  t>>A,  it  can  be  shown  that  the  third  term  of  (35) 
is  much  smaller  than  the  second  and  thus  can  be  ignored. 
Letting  A  approach  zero  (sampling  infinitely  fast),  (35) 
becomes 

P2  .  2W  )   sinfuWt  dt  (36) 

0    (irWt) 

and  then  since  t  is  large,  the  integral's  upper  limit  can 
be  replaced  by  infinity,  and  the  integral  reduces  to  one, 
giving  the  correct  result. 
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Actually,  by  including  the  third  term  as  A  goes  to  zero, 
P  will  be  slightly  less  than  one.   The  reason  for  this  is 
that  in  calculating  the  performance  figure  of  an  analog  TPR 
in  (7)  the  fact  that  -  <<  B  was  used  to  get  (5)  and  (6). 
Had  this  fact  not  been  used,  as  A  approached  zero,  F  would 
approach  one  exactly. 

P  will  be  plotted  as  a  function  of  WA.   Hence,  for  this 
case  the  Nyquist  sampling  rate  (twice  the  highest  frequency) 
corresponds  to  WA  =  1  . 

The  analysis  done  in  this  chapter  is  Just  a  special  case 
of  the  analysis  that  will  be  done  in  Chapter  V,  hence  all 
results  will  be  plotted  in  the  figures  of  that  chapter.   The 
results  obtained  from  (35)  for  an  ideal  bandpass  filter  are 
plotted  in  Figures  18  through  23.   It  is  of  interest  to  note 
that  for  W  less  than  one  (sampling  faster  than  the  Nyquist 
rate),  F  is  exactly  equal  to  one,  an  intuitively  satisfying 
result. 

2.   Single  Tuned  Filter 

For  a  Single  Tuned  Filter 


lHo(f)l2  = h^-;* b-r  "7) 

where  W  is  the  half  power  bandwidth  and  W  is  much  less  than 

f  .   It  follows  that 
o 

Sx(f)  =  -f   |Ho(f)|2  (38) 
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and 


and 


and 


N 

SY  (f)  =  Vt  (39) 

Xc      1  +  M§)2 

From  (39)  >  using  the  Inverse  Fourier  Transform 


N  W    ...  | 
Rx  (y)  =  -g—  e— n-Wlvl  (40) 


Top    T    T  k-1 


Equation  (.41)  can  be  put  in  closed  form  by  letting 
A  =  e"27rWA  (42) 


Using  the  fact  that  [Ref.  9] 


k=l       -1  "  A 


V  kAk  =  A-nAn+(n-l)An+1  (1|i() 

k=l  (l-AT 


and  after  some  tedius  but  straightforward  manipulations,  it 
follows  that 


2 

(fM2   =   f-  (coth  ttWA)    -  |p  (cosech2   irWA)(l   -  e"27rWx)      (45) 
op 

and 

F2   =   AB„(coth   ttWA)    -    (ABe)      fnnaa,   v,2   „taWi         q-2ttWTv    ,,.,-x 
E  id        (cosech     ttWA)(1   -  e  )    (^b; 

2TBE 
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It  can  be  shown  that  the  equivalent  bandwidth  of  a 
single  tuned  fulter  is  [Ref.  1] 

BE  =  ttW  (47) 

where  W  is  the  half-power  bandwidth.   Furthermore,  since 
t>>A,  the  second  term  is  negligible  with  respect  to  the 
first,  hence 

F2  =  7rwA(coth  ttWA)  (48) 

To  compare  with  the  analog  case,  let  A  go  to  zero. 
Then  P  would  become  one,  which  is  the  correct  result. 
Figures  24  through  29  show  the  results  obtained  for  this 
filter. 

3.   Gaussian  Filter 

For  a  Gaussian  Filter 

,     -v(£=£°)2     v(£l£2.)2 

[Ho(f)]2  =  e     W     +  e    W  (49) 

where  v  =  4  In  2  =  2.773,  W  is  the  half-power  bandwidth,  and 

W  is  much  less  than  f  . 

o 

It  follows  that 

-v(§)2 
SY..(f)  =  N  e   W  (50) 


and 


XcV  '    ° 


-!l!(Wy)2 
Rx  (y)  =  N0W(J-)J5e  v  (51) 
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Replacing  equation  (51)  Into  (30)  and  using  the  fact 
that  for  this  filter 

BE  =  (2JL)Sl  =   1.505  W  (52) 

it  follows  that 

2  n"1      kA    "|(lTkWA)2 

F^  =  1.505WA  +  3-01WA   E   (1-^F)  e  (53) 

k=l      T 

\ 
Using  the  same  argument  as  before  (A<<t),  (53)  becomes 

0  n-1    -|(TTkWA)2 

F^  =  1.505WA  +  3.01W  I   e  (54) 

k=l 


To  compare  to  the  analog  case  let  A  approach  zero.   Then 
(5^)  becomes 

P  t  -|(irkWA)2 

F^  =  3-01W  /  e  (55) 


and  since  x  is  large,  the  integral's  upper  limit  can  be 
replaced  by  infinity,  and  the  integral  reduces  to  one,  as 
it  should.   The  numerical  results  are  plotted  in  Figures  30 
through  35. 

.  4.   Butterworth  Filter 

The  last  filter  to  be  considered  is  a  Butterworth 
Filter  of  order  two.   For  this  case 

i  +  i6<-vr>      l  +  iec-jfV 

where  W  is  the  half-power  bandwidth,  and  W  is  much  less  than  f 
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Then 


N 

SY  (f)  =  2—r-E  (57) 

Xc      1  +  16(§)4 

and 

Rx  (y)  =  -V  e      (2)^  cos(l±^  _  J)  (58) 

Xc        2  (2)     4 


For  this  filter,  it  can  be  shown  that 

Bv   =  -^-^  ttW  =  1.48W  (59) 

E     3 

Substituting  (58)  and  (59)  into  (30)  it  follows, 
after  some  effort  that 


n— 1  h 

F2   =   1.48WA   +   2.96WA      Z   e~k7r(2)  2WA[1  +  sintku^^WA] 

k=l 

--2'96(AW)2     Vke^(2^WA[l  +  sin(k,(2)^A]      (60) 
TW  k=l 


Again  it  can  be  argued  that  the  third  term  is  negligible, 
therefore 


n— 1        Js 
F2  =  1.48WA  +  2.96WA   E   e'k7T(2)  2WA[1  +  sin(k7r (2)^2WA)  ]  (61) 

k=l 


Comparing  to  the  analog  case  by  letting  A  approach 
zero,  (6l)  becomes 

t      h 

F2  =  2.96WA  /  e"7r(2)2Wt[l  +  sin(TT(2)J'2Wt)]dt        (62) 
0 
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Since  x  is  large,  the  integral's  upper  limit  can  be  replaced 
by  infinity,  and  it  reduces  as  it  should  to  one. 

The  numerical  results  for  this  filter  are  plotted 
in  Figures  36  through  4l. 

B.   RC  FILTERING  AND  THEN  SAMPLING 

Refer  to  Figure  3.   This  type  of  digital  TPR  was  studied 
by  Ohlson  [Ref.  10]  for  the  case  where  no  quantizer  is 
present.   Only  the  results  and  its  implications  are  presented 
here. 

It  can  be  shown  that 

N 
E[Z]  =  g  -f  (63) 


and 


c72  =  g2  &)2   l (64) 


Z     &   v  2 


2BEtRC 


where  g  is  a  system  gain  parameter,  N  is  as  defined  on 
page  17,  and  tRC  is  the  time  constant  of  the  RC  filter. 
Also 


E[I]  =  E[Z]  (65) 


and 


crT2  =  -|-+2f   E   (1-  ^-)[R7(kA)-E2[Z]]    (66) 
1     n     T   k=l      T    L 


where  R7(u)  is  the  autocorrelation  function  of  the  process  Z(t). 
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Substituting  (65)  and  (66)  into  (13),  we  have 

_  ,  kA 

op       E  RC  k=l 

Defining 

A  =  e~A/tRC  (68) 

and  using  equations  (43)  and  (44)  it  can  be  shown  that 

\ 

(1      e-T/tRCw     A  A     s2 

P2    .       A       coth       A       _   ^J )^tRCCSCh2tRC)  (69) 

<^RC       ^RC  (T/tRC) 

Equation  (69)  is  of  the  same  form  as  equation  (47). 

Note  that  in  that  case  it  was  argued  that  the  second  term 

o 

was  negligible  since  tW  was  large  (oh  the  order  of  10  ) . 

Here  the  same  assumption  can  not  be  made  since  r — <<  W 

fcRC 

(i.e.  the  RC  filter  bandwidth  is  much  smaller  than  the 
bandwidth  of  the  RF,  Mixer,  IP  filter).  The  first  term  of 
equation  (69)  reduces  to  one  as  A  approaches  zero.   The 
second  term  however,  contributes  negatively.   Therefore  F 
can  be  less  than  one.   The  reason  for  this  lies  in  the  addi- 
tion of  an  RC  filter  to  the  radiometer ^  which  has  the  effect 
of  lengthening  the  integration  time  x.   As  A  approaches  zero, 
a  comparison  is  being  made  of  two  TPR,  one  with  an  RC  filter 
after  the  square  law  detector  and  the  other  without  one. 
The  one  with  the  RC  filter  will  have  greater  sensitivity 
(smaller  AT)  due  to  a  longer  effective  integration  time. 
Therefore  F  can  be  less  than  one. 
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However,  if  the  effects  of  sampling  only  are  of  interest 
(not  the  effect  of  the  addition  of  the  RC  filter),  the 
second  term  must  be  disregarded  by  assuming  that  T/tR„  =  °°. 
This  implies  that  all  the  filtering  is  done  by  the  integrator 
and  none  by  the  RC  filter. 

As  explained  in  part  A  of  this  chapter,  the  results  are 
presented  in  Chapter  V.   Refer  to  Figures  46  through  51. 

In  summary,  the  effects  of  an  ideal  ADC  (i.e.  no 
saturation  and  step  size  equal  to  zero)  was  studied  in  this 
chapter.   Its  main  importance  is  the  fact  that  it  provides 
a  lower  bound  for  the  Degradation  Factor  and,  that  under 
certain  conditions,  a  correction  can  be  added  to  these 
results  to  account  for  the  finite  step  size  of  the  quantizer. 
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IV.   TRANSFER  CHARACTERISTICS  OF  A  QUANTIZER 
WITH  A  ZERO  MEAN  GAUSSIAN  INPUT  PROCESS 

In  this  chapter  the  effects  of  a  quantizer  on  a  zero 
mean  Gaussian  input  process  are  analyzed.   Two  cases  are 
investigated.   Refer  to  Figure  6.   It  is  desired  to  find 
the  autocorrelation  function  at  point  2  given  the  auto- 
correlation function  of  X(t)  and  the  knowledge  that  it  is 
Gaussian  and  has  zero  mean.   Several  approaches  have  been 
explored  for  the  case  of  a  quantizer  alone.   The  one  sug- 
gested by  Kellog  [Ref.ll]  was  used  since  it  can  be  easily 
extended  for  the  case  of  a  quantizer  followed  by  a  square 
law  detector,  and  also  the  equations  that  it  leads  to  are 
In  a  form  suitable  for  computer  programming.   Figure  7 
shows  two  quantizers,  one  with  even  number  of  levels  and 
the  other  with  odd  number  of  levels. 

Let  N  be  the  number  of  levels  of  the  quantizer.   Then 


_  2L(i  -  N/2)  (?0) 

*i       n-1  UU; 


and 


n   =  2L[i  -  (N+l)/2]  rrn 

qi         N-1  w  J 
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A.       QUANTIZER   ALONE 

The   autocorrelation   function   of  q(t)    is 


Rq(y)    =   E{Q[X(t)]Q[X(t+y]} 


(72) 


Let 


X(t)    =    ?. 


(73) 


and 


X(t+y)    =   5, 


(74) 


Then 


00  00 


Rq(y)    =        f        f  Q(?1)Q(c2)f(?1,C2;y)dC1dC2  (75) 


_oo         _oo 
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where 


_(^-2pAC2-K22) 

(1-PX2)  2ax2 

f(C  .,?9;y)  =  - _ ^^  (76) 

1   2        a  2  2-rr  (l-pv2)* 


and 


RY(y) 
PX  "  •ixT0T-  (77) 


Noting  that  the  dependence  of  R  upon  u  is  only  via  pv, 

q  a 

the  notation  R„(PV)    will  be  used   for  brevity.      Then    (75) 
q       A 

becomes 


N        N  X.  X 

1-1   j-1  x±_1        Xj^ 


(78) 


Let 


B1  =  ?1/ax  (79) 


32  =  ?2/ax  (80) 


then    (78)    can  be  written  as 

h     i  -(B12-2pxe132+e22) 

ay       cty     ? 

N      N                                        p           2(1-0,-) 
Rq(px)  =     2      2  q1qJ  /  /       S . * d3id32        (8l) 

1=1  j=1  Xi-1      Vl  2(1  -p    2)^ 

_  _  A 

ax        ax 
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but 


Bl2    "   2f3ie2pX   +    322    =    (B2    "    PX31)2    +    Bl2(1   "    PX)2  (82) 


Therefore 

rh2   ~   2PXglg2   +    g22,  -    (g2    -    pXgl)2      612                   >Rq, 

-    L 5 J      =      5— 5—  (03; 

2(1    -   p/)  2(1    -    p/)           2 

and    (8lV  becomes 

p 

Xi  2                 X .       "(32"pX6l) 

N       N  „  n           ~ov       1  h  /2             oy    e  2(l-p  2) 

R  (py)  =     Z       E  qlqj       /     A    2~F    e  /    X X_d|3pd3..    (84) 

q    X        i-1  j=l          X^  X^      C2Tr(l-px2)]'2        2     X 

~°x~  aX 


Now  "let 


B2  -  p^ 

=  a 


^(l-px2)]"5 


(85) 


It    follows   that   the    second   Integral   can  be  written  as 


1  (erft  V"*  '  "ft]   -  erfC  i^X^  ^   ]}  (86) 

2  I  [2(1-Pv  )V  [2(1-PT2  )T  1 


where 


erf(a)    =   -4-       /      e"v      dv  (87) 

tt^      0 


35 


If  a  new   function   is   defined   as 

N    q.    r  X,/aY  -  py6,  X.   ,/o     -  p  8 

JL  l  ^ff  -J — £ £Ji  i   .Pr.f  r  -iz± — £ — JL±  ]    t     (88) 


X     X         j=l    2     l  [2(l-px2)]'2  [2(l-px2)] 

then    (84)    can  be  written   as 

Xi  2 

i       N  «F  -^ 

Rrt(py)    =     '        k        z      3,  /  ^(pY,en)    e      2        d(3n       (89) 

q     X  (2tt)^     i=l      1      Xi_1  X     X  1 

2 
Dividing  both  sides  of  (89)  by  ay  and  substituting  the 

expressions  for  X.  and  q. ,  it  is  found  that  a  new  parameter 

appears,  namely  aY/L  .   This  parameter  together  with  N  and 

Px  will  be  the  independent  variables  of  (89).   Note  that  with 

the  introduction  of  cry/L  ,  the  dynamic  range  of  the  quantizer, 

L,  is  no  longer  an  independent  factor. 

By  letting 

(90) 

(91) 


(92) 


qi  =   qi/L 

X£   =   X±/L 

and 

-    CTi  =  ax/L 

it   follows  that 

\ 

Vpx}  _     1 

^         (2ir)* 

n    q^ 
1=1  'X 

xi-i 

V 


^'(px,61)  e""2"  d61    (93) 


and 
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N  q»  f  X'/a'  -  py6n        X!  -,/a'  -  pi   ^ 

♦•(PjcA)  -  2  gX  erfC  J  X   M  ]  -  erf [  J"1  X   P  ]   (94) 


It  can  be  shown,  following  a  similar  procedure,  that  for 
p„  equal  to  one,  (93)  becomes 


Ra(l)  °a         1   N   qi  2  f     Xi/aX        Xi-l/aX   1 

-S_  =  J-  .  i   z  (  i. )2  erf[  -L-£  ]  -  erf[  -^-^  ]        (95) 

ax2    ax2   2  1=1  aX    I     (2)2         {2)h  > 

2 
This  gives  us  a   numerically  so  we  can  get  the  normalized 

autocorrelation  of  q.   Equations  (93)  through  (95)  allow  the 
computation  of  the  autocorrelation  function  of  the  output  of 
a  quantizer  of  N  steps  and  saturation  level  L  given  that  the 
input  is  Gaussian,  zero  mean,  and  of  known  autocorrelation 
function.   The  procedure  would  be  as  follows.   First, 
calculate  a  ' ,  X'  and  q'  for  the  given  quantizer  and  input 
process.   Second,  find  a  value  of  py  for  a  particular  u 
(Equation  (76)).   Third,  form  the  function  ^'(px,B1)  for 
the  particular  value  of  p„  determined  in  the  second  step. 
Fourth,  replace  the  function  ^'(px,B-.)  into  (93)  and  perform 
the  integration.   This  procedure  will  map  one  point  of  the 
input  autocorrelation  function  into  one  point  of  the  output 
autocorrelation  function. 

A  program  was  written  to  solve  (93)  through  (95)  for  the 
case  of  symmetrical  quantizers  (See  Appendix  under  Transfer 
Characteristics  of  Odd  level  Quantizer,  and  Transfer  Charac- 
teristics of  Even  level  Quantizer)  and  the  results  are 
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plotted  in  Figures  8  through  14  as  a  function  of  pv.   The 

Ra(v) 

ordinate  axis  is    a  0   =  p   .   In  other  words,  for  a  given 

0q2      q 

pY  ,   p   is  obtained.   The  advantage  of  plotting  it  this  way 
a     q 

is  that  the  curves  are  independent  of  the  input  autocorrela- 
tion function  (as  long  as  the  input  process  is  Gaussian  and 
zero  mean). 

In  general,  the  curves  are  a  function  of  tf-v/L  (except  for 

the  special  case  of  N=2).   Note  that  for  N=2  the  result  is 

\ 
the  well-known  relationship  [Ref.  7] 

MPX)    p 

q  0A  =  =-   arcsin  pY  (96) 

a   d  i\  A 

q 

It  is  of  interest  to  note  the  behavior  of  these  curves  as 
aY/L  approaches  extreme  values.   For  a  quantizer  with  odd 
number  of  levels,  as  tfy/L  gets  very  large,  the  quantizer  will 
appear  as  a  hard  limiter,  hence  the  curves  will  approach  that 
of  Figure  8.   As  ay/L  approaches  zero  the  quantizer  will 
appear  as  a  three  step  quantizer  (eventually  it  will  become 
an  open  circuit)  and  the  curves  will  show  a  behavior  similar 
to  those  shown  in  Figure  9.   For  a  quantizer  with  an  even 
number  of  levels,  the  behavior  of  the  curves  as  ay/L  becomes 
large  is  also  as  described  above.   When  cr  /L  approaches  zero, 
however,  the  quantizer  will  appear  again  as  a  hard  limiter, 
therefore  the  curves  will  approach  that  of  Figure  8. 
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B.   QUANTIZER  FOLLOWED  BY  A  SQUARE  LAW  DETECTOR 

The  reason  for  analyzing  this  case  will  become  apparent 
in  Chapter  V.   Refer  to  Figure  6  and  observe  that 

Y(t)  =  q2(t)  (97) 

and 

q(t)  =  Q[X(t)]  (98) 


R^y)  =  E[Y(t)Y(t+y)]  (99) 


Upon  substituting  (97)  and  (98)  into  (99)  it  follows 
that 

Ry(y)  =  E{Q2[X(t)]Q2[X(t+y)]}  (100) 

Following  a  similar  procedure  as  done  for  the  previous  case, 
(100)  becomes 

—  2 

Vfv)  i  N      <M  p        °x  "3l 

1    uA     =  — ±-r-       S      (rr)  /  A(pY,3n)    e~2~  dB,  (101) 

cx4  (2ir)*     i=l      ffX       Xjj  X     X  X 

~3T 


and 


N   q'     r         X'/a'  -  p  6n  X'  -,/a'  -  pY^ 


ACp-A  )  -  i  E  ftp  {erf  [    3     X    2  V  ]  -  erf  [  'J'1    X  "    P  ]}         ( 102 ) 
j-l°X  C2(1-PX2)]^  [2(l-px2)]^         J 

For   PY=1   it    can  ^e   shown  that 

ECY2]  _  1      ?     ,qix4  r     -  ,       Xi       ,  ...  Xj-i     v   -i  Mn,v 

-5-TT=?      E     Czt")     [erf  (  r  )  -  erf  ( j-  )   J  (103) 

V        d     1=1    aX  ax(2)^  ax(2)^ 
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Xi>  qi>  °X>    and  PX  were  deflned  ln  (90),  (91),  (92)  and  (77) 
respectively. 

A  computer  program  was  written  for  this  case  when  the 
quantizer  Is  symmetrical  (See  Appendix  under  "Transfer 
Characteristics  of  Odd  Quantizer  Followed  by  a  Square  Law 
Detector"  and  "Transfer  Characteristics  of  Even  Quantizer 
Followed  by  a  Square  Law  Detector").   The  same  comments  as 

in  the  previous  section  apply  for  the  plots  of  this  case 

RY(y) 
except  that  now  the  ordinate  axis  is   e — r^r  instead  of  py 

since  Y(t)  is  not  zero  mean.   The  results  are  shown  on 

Figure  15  for  N  equal  to  five  and  eleven. 

Ry(p) 

It  can  be  shown  that  for  the  case  of  no  quantizer  = — T-=-y 
is  given  by 

Rv(y)    1  +  2p„2 

^m-—r^-  (104) 

Note  how  the  curves  of  Figure  15  depart  from  the  curve  given 
by  (104)  due  to  the  effects  of  the  quantizer. 

In  summary,  in  this  chapter  an  algorithm  has  been  devel- 
oped suitable  for  computer  programming  in  order  to  find  the 
autocorrelation  function  of  the  output  of  a  quantizer  (or 
quantizer  and  a  square  law  detector)  given  that  the  input 
process  is  Gaussian  and  zero  mean,  and  with  known  auto- 
correlation function.   What  it  basically  does  (for  both  cases) 
is  to  reduce  a  double  integral  (75)  and  (100)  into  a  single 
integral  by  using  the  error  function  which  is  available  in 
most  computer  libraries. 
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V 


The  transfer  curves  for  both  cases,  shown  on  Figures  8 
through  14  and  15  are  independent  of  the  input  process 
providing  it  is  Gaussian  and  has  zero  mean.   These  results 
will  be  used  in  Chapter  V. 
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V.   DEGRADATION  FACTOR  CALCULATION  WITH 
BOTH  SAMPLING  AND  QUANTIZATION 

In  this1 chapter  the  effects  of  sampling  and  quantization 
on  the  performance  figure  of  the  two  digital  radiometers 
under  study  are  investigated.   As  in  Chapter  III,  the  two 
cases  will  be  considered  separately. 

A.   IF  SAMPLING 

1.   Exact  Results 


Refer  to  Figure  2.   In  Chapter  III  an  expression  was 

developed  for  the  degradation  factor  for  the  case  where  there 

was  no  quantizer  —  see  (30).   Those  results  cannot  be  applied 

directly  to  this  case  since  in  order  to  obtain  (30)  it  was 

assumed  that  q  (kA)  and  q  (kA)  were  samples  of  Gaussian 
c  s 

processes  —  see  (22).   With  the  addition  of  a  quantizer  it 
is  clear  that  this  assumption  no  longer  holds.   Hence  we 
must  start  again  with  a  more  general  approach.   This  can  be 
accomplished  by  concentrating  first  on  the  top  channel  of 
the  radiometer  of  Figure  2.   The  processing  that  the  ADC  and 
computer  do  to  the  signal  can  be  drawn  sequentially  as  in 
Figure  16.   It  is  argued  that  the  order  in  which  the  signal 
q  (t)  is  sampled  and  squared  is  immaterial  as  far  as  the 
random  variable  I,  is  concerned.   Hence  Figure  16  is  a  valid 
model  of  a  channel  of  the  radiometer  of  Figure  2. 

'  Referring  to  Figure  16  it  follows  that  [Ref.  12] 

a2        n-1 

a2  =  JL  +  I       i      (l-£)[RY(kA)  -  E2[Y]]  (105) 

±1    n      n  k=l     n   i 
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R„(k  )  is  the  k —  sample  of  the  autocorrelation  function  of 

Y(t) 

and 


E[I,]   =   E[Y(t)]   =   al 


(106) 


XcCO 


%l$\ 

LAW 
DETECTOR, 

Y(t) 

FIG! 

JRE   16 

L\Xc^|   ,^(^ 


Vc=» 


Using  (26)  and  the  fact  that  for  the  case  of  Figure 
2  the  RF,  Mixer,  IF  filter  is  symmetrical  about  f  ,  then 


«2  o    «2 

°1      =  2  °I. 


(107) 


and 


E[I]  =  2E[I,]  =  2  o^ 


(108) 


Applying  the  definition  of  degradation  factor  it 
follows  that 


2a: 


2   °Xc 


4BEA  n_1       kA  4 

o4   *"   <#   k=i      t    ^       q 


[dE£ll]2 
dot 


(109) 


aX    atTop 
c 
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Using  the  same  argument  as  in  Chapter  III  (A<<t)  it 
can  be  shown  that  the  term 


4BEA 


n-1 


kA 


2   ^     CRY(kA) 

k=l   T     1 


Mc 


(110) 


can  be  disregarded  due  to  its  small  value  compared  with  the 
other  terms  in  the  numerator  of  (110). 
Then  (109)  becomes 


P^  = 


2c%  4Bt.A  n-1  h 

ct  V+sdr    *     [RY(kA)  -  V 

Xr>       Xr«    k=l  ^C 


[  ffiEl  ]2 
daf: 


(111) 


at  T 


op 


Looking  at  (111)  term  by  term,  it  is  noted  that  in 

the  numerator,  the  first  term  can  be  obtained  from  (103)  and 

(95).   The  second  term  can  be  obtained  from  the  corresponding 

Rxc(kA) 

value  of  5 — /A\  .  and  the  transfer  characteristics  of  a 
RXC(°) 

quantizer  followed  by  a  square  law  detector  (Figure  15). 


Actually,  a  scale  factor  of  av  /ov   will  appear  because  of 

the  way  the  curves  were  normalized.   The  denominator  can  be 

obtained  from  (95).   Note  however  that  Equations  (93) »    (95) 

and  the  transfer  characteristics  developed  In  Chapter  IV  have 

aY/L  as  a  parameter  and  not  crv   .   This  is  due  to  the  choice 

A  Xc 

of  crx/L  as  a  parameter  since  it  was  found  in  that  chapter 

that  the  ratio  of  a     and  L  was  of  importance  and  not  their 
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absolute  magnitudes.   Therefore  it  makes  sense  to  have  (111) 

be  a  function  of  av  /L  at  T   instead  of  aY  .   av  /L  can 

Ae  op  Ac  Ac 

be  thought  of  as  a  normalized  standard  deviation  of  X  (t)  and 

its  value  can  be  related  to  temperature  in  the  same  form  as 

Equations  (10)  and  (11). 

2  2    2 

A  plot  of  E[I]/L  versus  cv  /L   is  shown  in 

Ac 

Figure  17  for  different  values  of  N.   The  effects  of  satura- 
tion of  the  quantizer  can  be  observed  as  temperature  (equiva- 

2    2 
lent  to  av  /  L  )  is  increased. 

Ac 

The  four  filters  previously  considered  for  this  case 

were  used  again  and  the  results  are  shown  in  Figures  18 

through  kl.      The  bottom  curve  for  each  figure  represents  the 

case  investigated  in  Chapter  III  where  there  is  no  quantizer. 

The  programs  used  to  evaluate  (111)  are  included  in  the 

Appendix.   It  should  be  pointed  out  that  the  programs  used 

to  calculate  the  transfer  characteristics  of  Chapter  IV  did 

not  provide  a  continuous  function  as  it  is  shown  on  the  plots 

of  that  chapter.   Since  in  order  to  solve  (111)  the  transfer 

curves  must  be  known  at  every  point,  a  least-mean-squared- 

error  polynomial  fit  was  done  to  the  data  obtained  from  the 

programs  of  Chapter  IV. 

Looking  at  Figures  18  through  Hi   it  is  noted  that  in 

each  case  there  is  a  reasonable  optimum  value  of  aY  /L 

Ac 

where  the  degradation  factor  is  lowest.   Physically  this  can 

be  explained  by  the  fact  that  for  small  a  /L   ,  the  step 

Ac 

size  degrades  the  performance  of  the  system,  hence  increasing 

aY  /L  will  tend  to  minimize  this  effect.   However  a  point  is 
Ac 
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reached  where  saturation  effects  start  dominating  and  the 

further  increase  of  a  /L  causes  an  increase  of  the 

degradation  factor.   For  N  >  4  ,   av  /L  optimum  is  around 

Ac 

0.35  regardless  of  the  sampling  rate.   This  would  imply  that 

for  a  given  T   ,  the  gain  setting  of  the  RF,  Mixer  IF  filter 

should  be  set  at  a  value  such  that   av  /L  is  equal  to  0.35 

Ac 

for  optimum  performance.   Note  that  this  optimality  is  for 
small  source  temperatures.   For  large  signals,  dynamic  range 
is  a  problem.   See  Chapter  VI. 
2.   An  Approximation 

The  evaluation  of  (111)  for  the  exact  value  of  the 
degradation  factor  requires  extensive  computer  time.   For  a 
quantizer  of  five  bits  or  greater  the  computation  time  per 
case  is  very  long  (on  the  order  of  hours).   Therefore  it  is 
of  interest  to  investigate  an  approximate  solution  for  the 
degradation  factor  and  to  compare  the  results  with  the  exact 
values  already  obtained. 

Figure  42  shows  one  of  the  channels  of  the  radiometer 
of  Figure  2.   Assume  that  the  sampling  rate  is  not  too  fast, 
the  step  size  of  the  quantizer  is  small  with  respect  to  the 
standard  deviation  of  the  process  X  (t)  and  no  saturation 

I 

takes  place  in  the  quantizer.   Then 

qc(kA)  =  Xc(kA)  +  ek  (112) 

where  X  (kA)  is  the  k —  sample  of  X  (t)  and  e.  Is  a  uniformly 
distributed  random  variable  that  accounts  for  the  error 
introduced  by  the  quantizer.   The  probability  density 
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function  of  ek  is  the  following  [Ref.  13] 


f(ek)    = 


1 
a 

0 


-  S.  <  e     <  £ 
2        ek        2 


(113) 


Then 


n 


n 


Xl  =  n"     kfx  qc^kA)    '   n-     k^    [Xc(kA)    +   2ekVkA)    +   ek] 


(114) 


and 


E[I1] 


i     Z    ECX^(kA)]    +  .£■      Z     E[2e.Xp(kA)]  +  i    Z     E[e£] 
n  k=]_        c  n    k=1  k  c  n  k=1        k 


(115) 


^Mf 


k-i 
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Assuming  e,  and  X_(kA)  are  uncorrelated,  which  is  a 

good  assumption  if  the  step  size  is  small  compared  with  o      . 

Xc 

2    a2 
E[I1]  =  ^x  +  12  (116) 

c 


Now 


E[I2]  =  -4-E  [  £  x2(jA)  +  2  E  e.X  (jA)  +  E  e2  ]  • 
n    j=l  C       j=l  J  C      j=l  J 


(117) 


[  E  XT(kA)  +  2  E  e.  X  (kA)    E  e2  ] 
k=l  c       k=l  K  c      k=l  * 


Expanding  (117)  and  using  (23)  and  the  assumption 

that  e,  and  e.  are  independent  for  k  /  j,  which  is  good  if 
K      J 

the  sampling  rate  is  slow,  it  follows  that 

"f,  =  ^r  +  4  j,  c**>4  (k»)  *jr  £  4.  +  J  A         <118> 


It  was  shown  previously  that  for  the  case  of  Figure 
2,  if  H  (f)  was  symmetrical  about  f  ,  then 

a2  =  2a2  (119) 

x     xl 

and 

ECI]  =  2  E[I1]  (120) 
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Using  (119)  and  (120)  it  follows  that 


2V  "Z1  „   k_2  „.,  ,  1  a2  „.   a4 


^  =  BEA  +  -?-   *  (1  "  ^  ^  CkA)  +  F  "IT  V  + -T-  350  V    (121) 


ax   ^="'"        c      0x      ay 
c  c       c 


But  the  first  two  terms  are  just  the  degradation 
factor  for  the  case  of  no  quantizer,  therefore  (121)  can  be 
written  as 

*2  -  PN«  *  \  4-   V  +  350  TT  V  <122> 

where  PNQ  is  the  degradation  factor  with  no  quantizer.   The 
last  two  terms  of  (122)  represent  a  correction,  due  to 
quantization,  for  the  degradation  factor  with  no  quantizer. 

The  approximate  solution  is  very  good  if  the  assump- 
tions under  which  it  was  derived  hold.   For  N=17,  o^   /L  =  0.3 


c 
(a/cr„  =  0.42)  and  WA=  1.0  the  error  is  less  than  one  percent. 
c 


For  values  of  a  /L  greater  than  0.3  the  approximation  is  not 

good  due  to  saturation.   For  values  much  smaller,  it  again 

fails  due  to  the  large  value  of  a/av  .   Figure  43  shows  how 

the  two  results  compare  for  N=7  and  different  values  of 

a„   /L.   In  it,  the  behavior  of  the  two  'solutions  as  o     /L   is 

c  !  c 

increased  can  be  observed.   For  N  much  greater  than  7,  the 

approximation  is  very  good  when  ov  /L  is  near  0.3. 

*c 
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B.   RC  FILTERING  THEN  SAMPLING 

Refer  to  Figure  3.   Unfortunately  there  is  no  exact 
solution  for  this  case  since  there  is  no  way  to  get  the 
second  order  probability  density  function  of  the  process 
z(t).   This  would  be  needed  to  approach  the  problem  in  the 
same  way  as  done  for  the  radiometer  of  Figure  2.   Two 
approximations  will  be  discussed. 

1.   First  Approximation 

It  was  shown  in  Chapter  III  that 


and 


E[Z]  =  g(NQ/2)  (123) 


4   -   S2<V2)2  2B^  <124> 


Owing  to  the  Central  Limit  Theorem  It  can  be  assumed 
that  z(t)  is  a  Gaussian  process  if  BptRr,  >>  1  .   Due  to  the 
square  law  detector,  it  is  obvious  that  z(t)  can  never  be 
negative,  hence  the  assumption  is  good  only  if  the  mean  of 
z(t)  is  sufficiently  large  compared  to  its  standard 
deviation  so  that  the  assumed  Gaussian  probability  density 
function  Is  negligible  for  negative  arguments. 

It  is  also  assumed  that  the  quantizer  has  a  bias  'b' 
as  shown  on  Figure  44.   When  b=0  this  quantizer  reduces  to 
the  one  shown  in  Figure  7.   It  is  argued  that  the  order  in 
which  z(t)  is  sampled  and  quantized  is  immaterial  as  far  as 
the  output  I  is  concerned.   Hence  Figure  45  is  a  valid  model 
of  the  processing  done  by  the  ADC  and  computer. 
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Let 


E[Z]   =  m 


(125) 


then  the  probability  density  function  of  z(t)  is 


fz(z)  = 


(2Tra?)?2 


(z-m)' 
o    2 


(126) 


It   can  be   shown,    using  a  procedure    similar  to  that 
introduced  in   Chapter  IV  that 


N      q,  X.H7} 

E[q(t)]  =     E      -±    [  erf(-    J 


i=l 


(2)% 


X.  --m 
.)  _  erf  (-±4— )   ] 


(127) 


Z 
Xj-m 


(2)% 


Vpz>  " 


N 
E     q. 


-32/2 


(2tt)^    i-1    X    Xj^-m  L 


dB 


(128) 


N      qj 

j=i  2  r 


iKp~,B)  =     E      4-    |erf[ 


X^-m 


-  P^B 


Vl"*1 


^(l-p2,)]1* 


]  -  erf [ 


0„-      ~pz6 


^(l-p2)]1* 


1 


]  (129) 
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and 


But  from  Figure  44 


X±  =  21.(1  -  |)/(H-1)  +  b  (130) 


q±  =  2L[i  -  (N+l)/2]/(N-l)  (131) 


Equation  (130)  can  be  written  as 


X,  =  X   +  b  (132) 

1    o± 


where  X   is  the  value  of  X.  when  the  quantizer  is  zero 

°i  i 

centered. 

Substituting  (131)  and  (132)  into  Equations  (127) 
through  (129),  the  parameter  b-m  appears.   Physically,  this 
is  the  offset  of  the  mean  of  the  process  z(t)  with  respect 
to  the  center  of  the  quantizer.   When  no  source  is  present 
the  value  of  this  offset  will  be  determined  by  the  quiescent 
conditions  of  the  radiometer.   When  looking  at  a  source, 
however,  the  value  of  the  offset  will  change  since  the  mean 
of  z(t)  changes  linearly  with  temperature. 

Let 

f  =  mQ  -  b  {133) 

where  m  is  the  expected  value  of  z(t)  at  T   .   Then  (127) 
through  (129)  can  be  written  as 
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N      qj   ,  X^  ~  f «  H2Btm)h  (az0  -  az) 


=  a7       Z      —  {  erf  [  r ] 

Z    i=l    a'    I  (2)^' 


E[q(t)]  _ 
L 

"'      ''Z    *  ^'  UZ 


(134) 


.  erf[    -i T    gj  "zo  -  ^  3} 


^i-l  "  f'  +(2BtRC^  (°Z0  "  az)  ^ 
(2)%a' 


'Z 


az' 


Rq(pZ)  j         N        q»  -B2/2 

-S5 h-    2        jr         ,      '  *,(pz,B)e  d6  (135) 

T2  (2ir)*i-l      aZ       Xo     i-f  Z 


*'(p7,B)  =    1      -1-    JerfC  -£ 5-t-  ]  -  erf  [  £ s-r—  ]       (136) 

Z  j=l    2az     I  C2(l-<72)]2  [2(l-a2)T*        J 


where  a „     is   the   value   of  the   standard  deviation   of  z(t)    at 
T      ,    and  all  primed  variables   imply  that   they  have  been 
normalized  by  L. 

For  the   case   of    p„  =   1      (135)    becomes 

E[-  2n        N    q'2  Xq.   -  f '  XL       -  f' 

EO  =     S      i       [erf  (  -2*-  )  -  erf  (       1-1  ^       )]  (137) 

c^  i=l  2cz  az(2)^  az(2)^ 


Equations    (13*0    through    (137)    can  be   solved  in  the 
same   form  as   done   in   Chapter   III.      However,   the   addition   of 
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two  new  parameters,  (f  and  (2BtRC)'5)  ,  makes  the  problem 
difficult  to  solve  in  a  general  way  due  to  the  number  of 
variables  involved. 

One  way  to  get  around  this  problem  is  to  assume  that 
(2B„tnr,P  is  large.   Then 


'E"RC- 


N        Xi 


E[q]  =   Z   q.     /   f7(z,m,a)  dz  (138) 

1=1      Xi-1 

where 

(z-m)2 
~   2 
f7(z,m,a)  =  K-r-     e   2a  (139) 


From  (123)  and  (124) 

m  =  ozC2BEtRC)'s  (140) 


Therefore 
-,„    N       X.    df7(z,m,a)  dz 

i-   *   11   /     -Lsa (m) 


but 
df 


Z  _  9f  +  3f  da  (142) 


dm    9m   9a  dm 

Substituting  into  (l4l)  and  using  the  fact  that 

55T  "  -  If  (l43) 

9m      9z 
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it  follows  that 


,  N  N     X  df 

f  =  ^W  *    qi[fz(Xi_1,m,a)-fz(Xi,m,a)]  +  Z    q±  Z1  ^  da  (144) 

X-  X  X- X     A.  -■ 


1   / 


% 


If  (2BFtRC)   Is  large,  it  means  that  the  expected 
value  of  z(t)  is  much  greater  in  value  than  its  standard 
deviation.   Therefore  it  can  be  assumed  that  a  change  in 
temperature  affects  the  expected  value  of  z(t)  but  the 
effects  onff„  are  negligible.   In  that  case  the  second  term 
of  (144)  can  be  disregarded. 

Using  (66)    it  follows  that 


°I:   ao    2A  n"1       kA   Ra(kA)  "  E  [q] 

4=^-+¥      E  u  -  ^H  - 5 1        0145) 

°Z  Z       k=1  °Z 

Substituting  (144)  and  (145)  into  (13),  it  is  easily 
shown  that 


a' 


«    ^RC      A    '  +  _±         A      ^        Rqto)  -  E2M 
.2  2        t™        ^ETRC  ~t™    ,  \     L  _2  J 


o      Gz  RC  RC    k=l  az 

(  dECq]  )2 
1    daZ     J 


2BE(A/tRC)2       n-1  Rq(kA)  -  E2[q] 

T/rRC  k=l  az 

,  dE[q]  v2 


(146) 
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The  parameter  2BEtRC  on  the  numerator  of  (146)  will 
cancel  out.   The  third  term  of  (146)  can  be  neglected  only 
if  T>>tR„  (see  comments  in  Chapter  III.   The  addition  of  a 
quantizer  does  not  invalidate  them). 

It  must  be  pointed  out  that  the  parameter  (2BptRC)  , 
even  though  it  does  not  appear  in  (146),  is  an  important  one. 
It  is  a  measure  of  how  much  smoothing  the  RC  filter  does. 
A  large  value  makes  quantization  difficult.   On  the  other 
hand,  a  small  value  invalidates  the  assumption  that  z(t)  is 
a  Gaussian  process.   Other  consequences  of  this  parameter 
will  be  discussed  in  Chapter  VI. 

Equation  (146)  can  be  evaluated  using  (13*0  through 
(137)  and  (144).   Basically  it  is  done  in  the  same  form  as 
for  the  case  discussed  in  Section  A  of  this  chapter.   Using 
(135)  and  (136)  a  set  of  transfer  characteristics  has  to  be 
computed  similar  to  Figures  8  through  14  (that  set  of  curves 
is  the  solution  when  f'=0).   The  major  difficulty  lies  in 
the  additional  parameter  f  which  increases  the  dimensionality 
of  the  problem. 

The  case  when  f'=0  (no  offset)  was  analyzed  using 
Figures  8  through  14  and  (144).   The  results  are  shown  in 
Figures  46  through  51.   For  any  other  value  of  f  (134) 
through  (137)  have  to  be  evaluated. 

As  an  example  the  case  of  N=13,  az/L=0.2  and  f'=0.5 
was  solved.   The  results  are  shown  in  Figures  52  and  53- 
Note  that  they  are  very  close  to  the  zero  offset  case  for  a 
seven  level  quantizer  with  az/L=0.4.   The  reason  is  that 
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a/az  is  the  same  for  both  cases.   They  are  not  exactly 
equal  however,  since  In  one  case  clipping  occurs  essentially 
for  positive  values  of  the  signal,  while  in  the  other  ■  case 
clipping  occurs  symmetrically  for  both  positive  and  negative 
values  of  signal.   Figure  52  is  a  plot  of  the  transfer  char- 
acteristics for  this  example.   The  programs  used  are  in  the 
appendix  (See  "Transfer  Characteristics  of  Quantizer  With 
Offset,"  and  "RC  Filtering  Then  Sampling,  Offset  Case"). 

The  case  when  N=2,  f'=0  is  of  interest  since  results 
are  obtained  easily  without  extensive  computer  programming. 
Using  (146)  and  (96)  it  can  be  shown  that 

o      w    a  a      n~1  VkA>         \r        n"1  VkA) 

^*4  +  4  A  "-jpbt  -it  &*"<** isr  (147) 

tRC 

The  last  term  can  be  neglected  when  x>>tRC.   As  the 
sampling  rate  increases,  A  goes  to  zero  and  Equation  (147) 
becomes 


F2  =  /  arcsin  e"y  dy  (148) 


For  large  x  ,  the  upper  limit  can  be  replaced  by 
infinity  and  (148)  integrates  to 

F2  =  \   In  2  (149) 


therefore 


F  =  1.04  (150) 
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The  results  are  surprising  since  it  means  that  a 
hard  limiter  degrades  the  system  only  by  four  percent. 
However,  the  dynamic  range  of  a  radiometer  using  a  hard 
limiter  would  be  very  small  hence  this  result  is  of  little 
practical  use.   More  on  the  problem  of  dynamic  range  will 
be  discussed  in  Chapter  VI. 

2.   Second  Approximation 

Using  the  same  arguments  as  in  part  A  of  this 
chapter,  it  can  be  shown  that  [Ref.  9] 

f2  -  pno  +  r — t~  (151) 

m        tRC      ct^24 

The  results  obtained  from  (151)  are  in  close  agree- 
ment with  those  obtained  using  the  first  approximation  even 
for  low  values  of  N.   It  must  be  kept  in  mind  that  the 
second  approximation  does  not  take  into  account  saturation 
effects  of-  the  quantizer.   For  the  case  when  f '=0  it  was 
found  that  if  a   /L  <  .6  WA  >  1  ,   a/a   <  2   and  N  >  5, 
the  second  approximation  was  within  1.5  percent  of  the  first 
(For  N=17  and  the  above  conditions,  the  difference  was  less 
than  0.6  percent).   These  results  indicate  that  a  sizable 
amount  of  saturation  can  occur  before  (151)  no  longer  holds. 
Figure  53  shows  how  the  two  approximation  compare  for  the 
example  done  on  page  91.   Figure  5^  shows  a  comparison  for 

N=7,  f'=0  and  different  values  of  a  /L.   The  different 

Z 

behavior  of  the  two  approximations  as  a„/L  is  increased  can 
be  observed. 
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To  summarize,  in  this  chapter  an  exact  solution  has 
been  found  for  degradation  factor  of  the  radiometer  of 
Figure  2.   For  the  radiometer  of  Figure  3  no  such  solution 
exists,  but  an  approximation  can  be  made  under  the  reasonably 
good  assumption  that  z(t)  is  a  Gaussian  random  process.   For 
both  cases  the  major  difficulty  in  global  analysis  lies  in 
the  number  of  parameters  involved. 

An  approximate  solution  was  developed  which  is  in 
close  agreement  with  the  first  solution.   Its  major  advantage 
is  that  it  is  easy  to  compute  and  as  the  number  of  steps 
increases,  it  gets  better.   For  practical  ADC  it  should  be 
adequate. 
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VI.   DYNAMIC  RANGE  AND  LINEARITY  OF  A  DIGITAL  RADIOMETER 

In  this  chapter  the  effects  of  the  addition  of  an  ADC 
on  the  linearity  of  theoutput  of  a  radiometer  is  investigated. 
The  main  problem  here  is  the  saturation  effect  of  the 
quantizer. 

A.   IF  SAMPLING 

Refer  to  Figure  2.   The  output  voltage  versus  temperature 
characteristics  for  this  case  was  investigated  previously, 
and  Is  shown  in  Figure  17.   If  the  ADC  was  ideal  (no  satura- 
tion and  a/a  =0)  there  would  be  one  straight-line  curve 
with  slope  equal  to  two.   It  can  be  seen  that  for  small 

av2/l2 (proportional  to  temperature)  the  curves  tend  to  have 
Ac 

a  slope  of  two. 

Figure  17  provides  all  the  information  needed  for  the 
problem  under  investigation  in  this  chapter.   However,  the 
way  the  curves  were  normalized  makes  them  awkward  to  use. 
They  can  be  changed  starting  with  (108)  and  (95). 

2 

-,rTn    aX    N  q'  2        X!  X'  ,  V^ 

^¥-   -  ~f   "  S  (=£-)  [erf  (-r-^-p  )  -  erf  (-r^-  1    (152) 

If  IT       1=1  aX  aX(2)  °X(2) 


but 


%       2 

TT  =  (N-l)a/ay  (153) 

Ac 
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Therefore    (152)   becomes 


E[I]       <4       N       q^(N-l)    2f  x;(N-l)  X^CN-1)        > 

-V"  =  — p-      s     [?CT     75  1   lerf[  r ]  -  erf[  p ]  f      (154) 

a2         a2     1=1    2aX/a        I         2(2)^av/a  2(2)^cv  /a     J 


—  Va  2u;ax 

c  c 


Equation    (154)    is   plotted  in  Figure   55.      The    curves   asymptote 
to 

sol  =  m£-  (155) 

a 

This  can  be  shown  by  using  the  fact  thas  as  cr„/a 
gets  large  the  quantizer  appears  as  a  hard  limiter. 

Figure  55  shows  that  a  substantial  range  of  linear 
output  can  be  obtained  with  this  scheme  before  saturation 
dominates. 

In  order  to  use  Figure  55  a  criterion  must  be 

defined  for  how  much  departure  from  linearity  is  allowed 

for  the  temperature  range  expected  (one  percent,  for  example). 

Once  that  is  defined,  Figure  55  can  be  entered  to  determine 

the  number  of  steps  of  the  quantizer  required.   The  value  of 

the  abscissa  at  T   will  be  determined  mainly  by  how  much 

op 

degradation  of  performance  is  allowed  (Chapter  V).   We  thus 
have  a  trade-off  of  linear  dynamic  range  versus  minimum 
detectable  temperature  AT. 
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B.   RC  FILTERING  THEN  SAMPLING 

Refer  to  Figure  3-   A  similar  analysis  can  be  made  for 
this  case.   However  the  additional  parameters  (f  and 
(2BEtRC)'1)  make  a  general  representation  difficult. 
Basically,  (127)  has  to  be  rewritten  in  terms  of  a„/a 
instead  of  a7/L.   It  can  be  shown  easily  that  (127)  becomes 


N  q!   ,    ^  -f'+(qgn^)s(2BEtRC)^  _ 

(2)%  s  a£ 


in  q:   , 
E[I]  =  o"       Z    -±-  4erf[ 
— T  1=1  2sa"  <• 


X£,   -  f  +  (ag  -<j")s(2BFt  )  2  ^ 

erf[  -^=1 j =2_Z ^RC   ]  i      (156) 

(2)^s  a£  _  J 


where  f ' ,  q!  and  X!  are  as  previously  defined,  and 

a£  =  az/a  (157) 


a"  =  a7/a  at  T  (158) 

Z0    Z         op 


s  =  2/(N-l) 


Figure  56  is  a  plot  of  (156)  for  the  case  of  f'=0, 

w 
a£  =0.5  and  (2Bc,tTDr,)   =  1.0.   The  abscissa  axis  is  shown 

as  a  function  of  crg/a^   (or  T/T   ) .   The  curves  asymptote 

o 
to  (N-D/2.   Some  observations  are  of  interest.   First,  the 

h 
slope  near  the  origin  is  equal  to  (2BEtRC)  a£   .   Since  for 

a  given  quantizer  E[I]/a  reaches  a  fixed  value  for  large 
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values  of  crj/ag   ,  increasing  (2BEtRC)   (which  physically 

o 
means  that  more  smoothing  will  be  done  by  the  RC  filter) 

decreases  the  range  of  linear  output.   Therefore  tR~  should 

be  chosen  as  the  minimum  value  which  will  enable  the  sampler 

to  operate  properly.   Second,  the  parameter  a1!     will  be 

Lo 

determined  by  the  amount  of  degradation  tolerated  at  T 

op 

And  third,. the  desired  range  of  linear  output  will  determine 
the  number  of  steps  of  the  quantizer.   Therefore  a  procedure 
to  use  a  Figure  such  as  56  would  be  as  follows:   First, 
determine  (2BEtRC)^  and  a"  using  the  criteria  given  above. 
Second,  draw  a  set  of  curves  of  the  expected  value  of  the 
output  versus  temperature  using  (156).   Third,  decide  upon 
a  criterion  of  how  much  departure  from  linearity  is  allowed 
in  the  temperature  range  of  operation,  and  finally  obtain 
the  required  number  of  steps  from  the  curves  drawn. 

To  summarize,  the  effects  of  the  addition  of  an  ADC  on 
the  linearity  of  the  output  of  the  radiometers  under 
consideration,  have  been  investigated  in  this  chapter.   The 
knowledge  of  the  behavior  of  the  output  is  important  for 
calibration  purposes. 
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calibration  purposes. 
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where  F  is  the  degradation  factor  for  the  equivalent  total 
power  radiometer.   The  factor  multiplying  F  is  the  performance 
figure  of  a  balanced  Dicke  radiometer. 


()-®-nVM»!- 


k--» 


FIGURE  57 

Equation  (160)  follows  from  the  fact  that  for  the  half 
of  the  switching  cycle  that  the  radiometer  is  connected  to 
the  source^  it  behaves  as  a  total  power  radiometer. 

Care  must  be  taken  however  on  the  selection  of  t.   If 
commencement  of  sampling  is  delayed  after  switching,  it  must 
be  corrected  for  the  loss  in  integration  time.   Also,  for 
the  case  of  RC  filtering  followed  by  sampling,  the  introduc- 
tion of  the  RC  filter  must  be  reflected  in  an  increase  of  t 
since  F  takes  only  into  account  the  effects  of  the  inclusion 
of  the  ADC. 

The  problem  of  dynamic  range  and  linearity  of  the  output 
can  be  handled  with  the  curves  developed  in  Chapter  VI. 
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B.   NOISE  ADDING  RADIOMETER 

Figures  59  and  60  show  two  digital  noise  adding  radio- 
meters.  The  one  in  Figure  59  corresponds  to  the  IF  sampling 
case  and  only  one  channel  is  shown.   The  one  in  Figure  60 
corresponds  to  the  RC  filtering  case.   The  extension  of 
results  already  found  to  NAR  is  considerably  more  difficult 
than  for  a  balanced  Dicke  radiometer.   This  is  due  to  the 
processing  involved  [Ref.  13] • 

Let 


and 


I  =  E  +  e  (161) 

n    n    n 


Io  =  EQ  +  eo  (162) 


where  I  and  I  are  the  outputs  of  the  radiometer  with 
no 

T  +T  and  T  _  respectively.   E_  and  E^  are  the  expected 
op  n      op  n      o 

value  of  the  output  for  the  case  of  T  +T  and  T   respec- 

*  op  n      op    * 

tively.   e  and  e  are  independent,  zero  mean  random  variables 

with  standard  deviation  equal  to  a,   and  aT   ,  and  represent 

n       o 
the  deviation  from  the  mean  for  a  particular  measurement 

(one  cycle). 

Then  the  Y-factor  is  [Ref.  14] 

en 

I    E^    1  +  p2- 

Y  =  -a  =  gH  [     En  ]  (163) 

o    o       e0 

1  +  ^ 
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If  en  and  eo  are  small  compared  with  E  and  E  ,  then 


E         E 

Y  =  —  [  l  +  Jl 

E^  L  x        E 


n 


E. 


(164) 


E 
E[Y]  =  =£ 
Eo 


(165) 


and 


a2  =  c  *a  ]2  C(in)2^,-(i°)2] 


n 


Using  Equation  (166)  and  the  fact  that  for  a  NAR 


(166) 


.AT.  .  (1  +  ^  <  _i_  ,* 

op         n      E 


(167) 


where  t  is  the  switching  period;   for  the  case  of  Figure  2 
it  can  be  shown  by  straight  substitution  that, 


1  r  ,  dE"n]  g 

dox 


NAR 


aX         dE[I  ]   2 

(1+^Fn  +  <  —2^-  >2 

X  X 

T   +T      °    ....  E[I  J. 
op  n o 


O 


^p_ 


(1  + 


2      d^n] 

o,2    r  .2 
— )   [  dcx 


n 


E[In] 


dE[Io] 


da' 


T  +T 
n  op 


etti 


n' 


"op 


(168) 


where  c\r  is  a,r  at  T  .,  ov  is  av   at  T 


X 


op'   X 


n 


X     n. 
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F  is  the  degradation  factor  of  the  corresponding  total 

power  radiometer  at  T=T  +T   ,  and  F  is  the  degradation 

n  op       o  ° 

factor  of  the  corresponding  total  power  radiometer  at   T=T 

op 

For  the  case  of  Figure  3,   o„   must  be  replaced  by  a„   in  (168). 

All  the  terms  of  (168)  can  be  calculated  from  equations 

previously  developed. 

It  must  be  noted  that  for  the  case  when  the  signal  is 

filtered  before  sampling  (Figure  60),  f  (offset)  cannot  be 

zero  at  T  ^  or  T  +T  since  the  ADC  would  give  zero  output, 
op     op  n  to  ^ 

The  problem  of  dynamic  range  and  linearity  is  easy  to 
handle.   From  Figures  55  or  56,  depending  on  the  system  under 
consideration,  and  (165)  a  plot  can  be  made  relating  E[I] 
to  temperature.   Figure  6l  shows  how  this  is  done  for  the 
system  of  Figure  59. 

In  summary,  in  this  chapter  the  results  obtained  in  the 
six  previous  chapters  for  the  digital  total  power  radiometers 
shown  in  Figures  2  and  3,  have  been  used  to  determine 
performance  of  other  digital  radiometers.   The  balanced 
Dicke  and  noise  adding  radiometers  were  investigated.   For 
the  balanced  Dicke  radiometer,  there  is  a  simple  and  straight- 
forward relationship.   For  the  NAR,  however,  a  relationship 
was  found  but  it  is  more  involved.   The  problem  of  dynamic 
range  and  linearity  can  be  handled  easily. 
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VIII.   DIGITAL  FILTERING 

In  the  analysis  of  the  classical  total  power  radiometer 
an  Integrator  is  used  to  smooth  the  output  signal.   This 
scheme  works  well  but  it  has  the  disadvantage  that  it  gives 
output  information  only  at  discrete  intervals  of  time,  corre- 
sponding to  blocks  of  data.   Such  performance  might  not  be 
acceptable  in  some  cases  because  a  smooth  output  might  be 
desirable.   Therefore,  a  pure  integrator  is  not  generally 
used  but  rather  a  low  pass  filter  (such  as  a  RC)  that  pro- 
vides a  continuous  output.   The  theory  and  equations  of  a 
classical  total  power  radiometer  do  not  change  when  using 
these  low  pass  filters  instead  of  an  integrator  if  an 
equivalent  integration  time  is  defined  [Ref.  $]. 

The  digital  total  power  radiometers  of  Figures  2  and  3 
suffer  the  same  disadvantage  as  the  classical  TPR.   That  is, 
output  information  is  obtained  only  after  the  summers  have 
completed  the  summation  of  the  n  samples.   Therefore  a  scheme 
where  output  information  is  obtained  more  often  must  be 
investigated.   This  chapter  deals  with  the  substitution  of 
the  summers  of  Figures  2  and  3  by  digital  filters.   There  is 
an  additional  reason  for  wanting  to  do  this  analysis.   Since 
a  computer  is  already  being  used  in  the  radiometers  of 
Figures  2  and  3,  the  use  of  a  digital  filter  often  does  not 
require  any  more  hardware. 

Figure  62  shows  a  digital  filter  operating  on  a  digital 
signal  X.  .   The  filter  is  characterized  by  its  impulse 
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FIGURE  62 


response  h,  ,  and  the  output  is  Y,  .   The  filter  is  assumed 
to  be  causal  (i.e.  h,  =0  for  k<0)  and  time  invariant.   This 
filter  can  replace  the  summers  of  the  radiometers  defined  in 
the  previous  chapters. 

Referring  then  to  Figure  62 


Yk  -  E  xi\-i 

i=-oo 


(169) 


By  straightforward  algebra  it  can  be  shown  that 


E[Y,  ]  =  E[Y]  =  E[X]   E   h, 
K  i=0  1 


(170) 


and 


o-J  =   Eh,   E  [RY    -  E^CXirjh^ 
Y   i=0   1  m=0  Xm-i         m 


(171) 


where  Rv  =  E[X,  X.  ,  ]  is  the  discrete  autocorrelation  of  X,  . 
X       k  k+n  k 

n 


Let 


J  =  i  -  m 


(172) 
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then 


Oy   =   E  h    E   [R   -  E*[X]]h,  ,  (173) 


Assume  that  h  decays  much  slower  than  Rv  .   This  can  be 
•J  Xj 

seen  by  observing  the  location  where  this  filter  would  go 
in  Figures  2  and  3-   For  the  case  of  Figure  2,  the  power 
spectral  density  at  the  input  of  the  digital  filters  is 
shaped  by  the  RF,  Mixer,  IF  filter  which  is  broad-band. 
For  the  case  of  Figure  3,  it  can  be  assumed  only  when  tRf, 
is  small  (this  is  the  case  for  which  Figures  46  through  51 
apply). 

Then  (173)  becomes 


oo  4 

2 
aY  =  ±lQ   hi   jf_„  CRXj  -  E2[X]]h±  (174) 


Let  j   be  the  value  where  R(j  )  is  negligible.   Then 
(17^)  can  be  rewritten  as 

o      i  °°       i 

E   h2    E   [R„  -  E2[X3]  +   E   h2    E  [RY  -  E2[X]] 
i=0   x   j—   Aj  i=i   1   1=-oo  A1 


a;  =  2  h 


(175) 


It  can  be  shown  that  the  first  summation  is  much  smaller 
than  the  second,  and  hence  it  can  be  ignored.   Furthermore 

the  limits  of  the  second  summation  can  be  changed  as  follows 

2 
without  altering  the  value  of  ay  significantly 
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a2  =   E  h2   E   [R   -  E2[X]] 
i=0     j=-°°   j 


(176) 


For  the  summer  used  in  previous  chapters,  it  can  be  shown 
(using  the  same  notation  of  Figure  62)  that 


a2  =  a2  +  2n   E   [R   -  E2[X]]-  2n   E   £  [Ry   -  E2[X]]   (177) 
1         A  k=l  Ak  k=l  n   *k 


and 


E[Y]  =  E[X] 


(178) 


The   upper   limit   of  the    summation  of   (177)    can  be   replaced 
by   infinity    (due   to  the   narrowness   of  RY,(u))    and  the   third 
term  can  be   ignored    (see   Chapters   III   and  V). 

Substituting   (170)    and    (176)    into    (4)    for  the    general 
digital   filter,    it    follows   that 


oF>2 

op 


o  °°  °° 

P       E    h.     +    p-      E     h. 

Tj_       i=0     J  T__      i=0     x    j=l    *j 


E  [Ry    - 


E2[X]] 


op 


op 


(zh)2     (dED0)2 
1-0 *  dT 


op 


(179) 


Doing  the    same  with    (177)    and    (178)    it    follows   that 


op 


n  T™2 


E     [Ry    - 


nV    k=l      *k 


E2[X]] 


rdECXl-,2 

L  dT     J 


(180) 


T 
op 
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Now  if  a  comparison  is  made  of  (179)  and  (180)  it  is 
clear  that  they  are  equal  if  np  is  defined  as  follows 


c   2  h,  )2 

-E  -   ■1:°  -  (18D 

Z     tu 
i=0   x 


and  n„  can  be  interpreted  as  an  equivalent  summation  number. 
Physically  it  is  a  measure  of  the  performance  of  a  general 
digital  filter  in  terms  of  the  performance  of  a  pure  summer. 
In  other  words,  it  provides  a  way  to  replace  the  digital 
filter  by  a  summer  of  np  samples  without  changing  the  per- 
formance of  the  radiometer.   Then  for  a  given  sampling  rate 
and  known  nF  an  equivalent  integration  time  can  be  determined 
and  a  performance  figure  can  be  calculated  using  the  curves 
of  Chapter  V. 

For  a  digital  total  power  radiometer  such  as  those  of 

AT 
Figures  2  or  3  using  a  digital  filter,  ss —  would  be  given  by 

op 


AT   _    1  (l82) 


Top    (BgngA)1* 


where  F  is  the  degradation  factor  for  the  equivalent 
radiometer  using  summers  and  (BEn„A)  2  is  the  performance 
figure  of  a  total  power  radiometer  with  Integration  time 
equal  to  n„A . 
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It  is  often  convenient  to  express  (l8l)  in  its  equivalent 
z-transform  since  the  transfer  function  of  a  digital  filter 
is  usually  given  in  terms  of  its  z-transform.   It  can  be 
shown  that  [Ref.  15] 


(H(z) 


nE  = 


)2 
2Zl (183) 


-In  -1 


-jp=-  <j>  H(z)H(z~x)z~xdz 

r 

where  the  contour  integration  denoted  by  r  is  around  a  unit 
circle  and  H(z)  is  the  transfer  function  of  the  filter. 

.'.  As  an  example  of  the  use  of  (183)  a  first  and  second 
order  digital  filter  will  be  considered. 

For  a  first  order  filter,  the  transfer  function  is 
given  by  [Ref.  13] 


H(z)  =  , ^__  (184) 

1  -  rz 


It  can  be  easily  shown  that  (183)  becomes 


by 


°E  -  TH  <l8*> 


The  transfer  function  of  a  second  order  filter  is  given 


H(z)  =  ZI± ^ (186) 

1  -  z   a  -  z   X 


Substituting  (186)  into  (183)  and  after  considerable  algebra, 
it  can  be  shown  that 
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„   _  X4  -  X2(g2+2)  -  2Xa2  -  a2  +  1  ,_«_. 

n„  -  5 5 (lo7) 

E     (1-a-X)2  (1-X2) 


Equation  (182)  and  (183)  are  very  similar  to  the  ones 
defining  the  equivalent  integration  time  of  a  low  pass  filter 
of  an  analog  total  power  radiometer.   This  is  not  surprising 
since  many  similarities  exist  between  digital  and  analog 
networks  even  though  they  are  not  equivalent. 
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IX.   CONCLUSION 

This  work  investigates  the  performance  of  radiometers 
using  digital  processing.   It  is  assumed  that  the  changing 
of  the  signal  from  an  analog  to  digital  form  is  performed 
before  the  final  output  is  obtained.   Two  digital  total 
power  radiometers  were  considered  in  detail. 

In  Chapter  II  a  figure  of  merit  is  defined  in  order  to 
evaluate  performance.   The  effects  of  the  addition  of  an 
analog  to  digital  converter  are  discussed. 

In  Chapter  III  the  effects  of  an  ideal  ADC  are  investi- 
gated.  It  serves  as  an  introduction  to  the  problem  of 
sampling  and  quantization.   The  results  are  useful  since 
they  provide  a  lower  bound  for  the  degradation  factor  as 
defined  in  Chapter  II. 

Chapter  IV  deals  exclusively  with  the  effects  of  a 
quantizer  on  the  autocorrelation  function  of  a  zero  mean 
Gaussian  input  process.   The  problem  of  a  quantizer  followed 
by  a  square  law  detector  is  investigated. 

Chapter  V  combines  the  information  of  Chapters  II,  III 
and  IV  and  treats  the  effects  of  a  real  ADC.   Both  quantiza- 
tion and  sampling  effects  are  investigated.   For  one  of  the 
radiometers  under  investigation,  exact  results  are  obtained. 
For  the  other  a  reasonably  good  assumption  is  made  and  an 
approximate  solution  is  obtained.   Other  approximation  are 
discussed  and  compared. 
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Chapter  VI  discusses  the  problem  of  dynamic  range  and 
linearity  involved  in  a  digital  radiometer. 

In  Chapter  VII  the  results  obtained  in  the  previous 
chapters  are  used  to  investigate  other  types  of  digital 
radiometers.   Balanced  Dicke  and  noise  adding  radiometers 
are  considered. 

Chapter  VIII  investigates  the  problem  of  using  a  general 
digital  filter  for  processing  of  signals  in  a  radiometer. 

The  major  difficulty  encountered  throughout  the  work  was 
the  number  of  variables  that  enter  into  the  problem.  There- 
fore it  was  Impossible  to  do  a  global  analysis.  However, 
the  factors  and  considerations  that  enter  into  the  design  of 
digital  radiometers  were  treated  in  sufficient  depth  so  that 
a  good  understanding  of  the  advantages,  problems  and 
limitations  of  digital  radiometers  Is  presented. 
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APPENDIX 

The  following  Is  a  list  of  the  programs  used  to 
evaluate  the  equations  and  draw  the  figures  throughout  the 
thesis.   They  are  listed  In  the  order  they  were  used. 


C      TRANSFER  CHARACTERISTICS  OF  ODD  LEVEL  QUANTIZER 

C      N=NUMBER  OF  STEPS 

C      VAR=NORMALIZED  STANDARD  DEVIATION  OF  X 

C      RHC=COPRELATION  COEFICIENT  OF  X 

C      DEL=INTEGRATION  STEP 

C      XF=X(N) 

C      XC=X(0) 

IMPLICIT  REAL*8(  A-H,C-Z) 

DIMENSION  SYS{ 5001) ,Q(34) ,  X(34) , FUN (5001) ,W(18,2) 
1  »  M  ( 1 8  .  2  ) 
C      PROBLEM  PARAMETERS 

READ(5,2C7)N 
207  F0RMATII2) 

WRITE(6,666)N 
666  FORMAT  ('  ',12,//) 

DO  205  NN=1,8 

READ(5.2G0)VAR 
200  F0RMAT<F4.1) 
C      ALGORITHM  PARAMETERS 

CEL=„001 

RH0=0.0 

R00=0.0 

XF=1000. 

XO=-1000. 

SY=-5.0 

SS  =  0.0 

S=0.0 
.  E=0.0 

STP=(2./(N-1.) ) /VAR  V 

PP=. 3989423 

M  =  (N-1  )/2 

L=DABS( SYJ/DEL 
C      QUANTIZER   CHARACTERISTICS 

DC  81  1=1, N 

0(I)=( I-CN+1.) /2.0)*STP 

X(  I)  =  { I-N/2.01*STP 
81  CONTINUE 
C      LIMITS  OF  INTEGRATION 

CO  97  J=1,L 

SYS{ J)=SY 

SY=SY+DEL 
97    CONTINUE 

DO    91    IR=1,N1 

IF(IR.NE.l)    GO    TO    92 

M(l,l )=1 

K3  =  (X(  D+5.0)  /DEL+1. 

WG=X( 1)-SYS(K3) 

M(  1,2)  =  K3 
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W(1,2)=W0 

GO    TO    91 
92    M(IR.1)=M(IR-1,2)+1 

W( IR,l)=DEL-WO 

K4=(X(  IRJ  +  5.0)  /DEL+1. 

WG=Xl IR )-SYS(K4) 

V  (  I  R  ,  2  )  =  K4 

W(  IR,2)  =  W0 
91     CONTINUE 

FINAL    POINT 

CO    98     1  =  1,  N 

IFU.NE.l)     GO    TO    99 

E2=X0/1.414 

E1=X( I J/1.414 

GO    TO    100 
99    IF( I.NE.N)     GO    TO    101 

E2  =  E1 

E1=XF/1.414 

GO    TO    100 

101  E2=E1 

E1=X(I J/1.414 
100    R2=QERF(E1)-DERF(E2) 
RR=Q( I )**2*R2 
SS='SS+RR 
98    CONTINUE 

RATT=SS/2.0 

RAT=DSORT(RATT) 

WRITE(6,102)RATT,"RAT 

102  F0PMAT(2F15.9) 
SY=-5.D 

DO    30    MM=1.21 

FORM    FUNCTION    PHI 

A  =  2.0*{  1.0-RHO**2) 

AA=DSORT(A) 

DO    2    J  =  1,L 

DO    1    I  =  1,N 

IF(I.NE.l)     GO    TO    5 

E2=(X0-RH0*SY) /AA 

E1=(X( I )-RHO*SY)/AA 

GO    TO    25 
5    IF( I.NE.N)     GO    TO    3 

E2=E1 

E1=(XF-RF0*SYJ/AA 

GO    TO    25 
3    E2  =  E1 

E1=(X(I )-RHC*SY)/AA 
25    R1=DERF(E1)-DERF(E2) 

R=0( I}*R1 

S  =  S+R 

1  CONTINUE 
FUN(J)=(S/2.0)*DEXP<[SY**2)/{-2.0)) 

.   S=.0 
SY=SY+DEL 

2  CONTINUE 

MAIN    INTEGRATION 

DO    7    11=1, Nl 
44    IF( II .NE.l)    GO    TO    11 

M1=M( 1,1) 

M2=M(1,2) 

C0RR  =  H(1  ,2)*FUN{M2) 

GO    TO    15 
11    M1=M(I1,1) 

M2=M(I1,2) 

CORR  =  W{  I  l,l)*FUN(Ml)  +  W(Il,2**FUNt-M2> 
15    M3=M2-1 

DO    21    K0=M1,M3 

011=FUN(KC)+FUN(K0+1) 

B=B+011 
21    CONTINUE 

R0=0( Il)*( DEL*B+CORR) 

FCQ=RQO+RO 

B=.0 
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7    CONTINUE 
PCOO=ROC*PP 
PGQO=(l. 3/RATT)*RQGC 
WRITE(6,40)RHQ,ROQQtRQQO 
40    FORMAT  (3D25. 9) 

WRITE( 7,204) RQ 00, ROCG 

204  FORMAT (2F8. 5) 
RCQ=0.0 

s=.o 

SY=-5.0 

IFtMM.EQ.20)    GO    TO    67 

RH0=RHO+.05 

GO    TO    30 
67     PH0=RH0+.03 
30    CONTINUE 

205  CONTINUE 
STOP 

END 


C 

c 
c 
c 
c 
c 
c 


666 


200 


81 


97 


92 


TRANSF 
N=NUMB 
VAR=NO 
RHC=CO 
DEL=IN 
XF=X(N 
XO=X(0 
IMPLIC 
DIMENS 
1,M(18, 
PROBLE 
N=4 
WRITE! 
FORMAT 
DO  205 
READI5 
FORMAT 
ALGORI 
RHO=D. 
CEL=0. 
RCQ=0. 
XF=100 
XC=-10 
SY=-5. 
SS=0.0 
S=0.0 
B  =  0.0 
STP=(2 
PP-.39 
Nl=N/2 
L=CABS 
OUANTI 
CO  81 
0(I)=( 
X(I)=( 
CCNTIN 
LIMITS 
DO  97 
SYS( J) 
SY=SY+ 
CONTIN 
DO  91 
IF1IR. 
M(l,l) 
K3=(X( 
WO=X(  1 
M( 1,2) 
W(l,2) 
GO  TO 
M{  IR.l 


ER 

ER 

RM 

RR 

TE 

) 

) 

IT 

10 

2) 

M 


CHARACTERISTICS  OF  EVEN  LEVEL  QUANTIZER 

OF  STEPS 
ALIZED  STANDARD  DEVIATION  OF  X 
ELATION  COEFICIENT  OF  X 
GRATIOM  STEP 


REAL«8(A-H,0-Z) 
N  SYS( 5001) ,0(34) ,X( 34)  , FUN  I  5001)  ,W(18,2) 

PARAMETERS 


6,666) M 

('  «,I2,//) 

NN=1,8 
,200) VAR 
(F4.1) 

THM  PARAMETERS 
0 

001 
0 
0. 
00. 
0 


./(N-l.)  JVVAR 
89423 

(SYJ/DEL+1 

ZER   CHARACTERISTICS 

1  =  1, N 

I-lN  +  1  .)/2.0)*STP 

I-N/2.0)^STP 

UE 

GF  INTEGRATION 
J  =  1,L 
=  SY 
DEL 
UE 

IR=1,N1 

NE.l)  GO  TO  92 
=  1 

1 )+5.0)/DEL+l. 
)-SYS(K3 ) 
=  K3 
=  WC 
91 
)=M( IR-1,2)+1 
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W( IR.l )=DEL-WO 

K4=(XUR)  +  5.0)/DEL  +  l  . 

WO=X(IR )-SYS(K4) 

M( IR,2)=K4 

W(  IR  .2)=W0 
91  CONTINUE 

FINAL  POINT 

DO  98  1=1. N 

IF(I.NE.l)  GO  TO  99 

E2=XC/1.414 

E1=X{  I  )  /1.414 

GO  TO  100  . 
99  IFfl.NE.N)  GC  TO  101 

E2=E1 

E1=XF/1.414 

GO  TO  100 

101  E2=E1 

F1=X( I  1/1.414 
100  R2=DERF( E1)-DERF(E2) 
RR=Q( I )**2*R2 
SS=SS+PR 
98  CONTINUE 

RATT=SS/2.0 

RAT=DSQRT(RATT) 

WRITE{6,102)RATT,RAT 

102  F0RMATC2F15.9) 
SY=-5.0 

DO    30    MM=1,21 
FORM    FUNCTION    PHI 
A=2.0*(1.0-RHO**21 
AA=DSORT(A) 
CO    2    J=1,L 
DO    1    1=1, N 
IFCI.NE.l)    GO    TO    5 
E2 = (XO-R HO* SY)/AA 
El=( X( I )-RHC*SY)/AA 
GO    TO    25 
5     IF  ( I  .NE.M     GO    TO    3 
E2  =  E1 

E1=(XF-RH0*SY) /AA 
GC    TO    25 
3    E2=E1 

E1=(X( I )-RHO*SY)/AA 
25    R1  =  DERF(E1)-DERF(E2  ) 
R=0(I )*R1 
S=S+R 

1  CCNTINUE 

FUN(J)=( S/2.0)*DEXP(( SY**2 ) / (-2. 0) ) 

S=.0 

SY=SY+OEL 

2  CONTINUE 

•   MAIN    INTEGRATION 

DO    7    I1=1,N1 
44    IF(Il.NE.l)    GO    TO    11 

M1=M(1,1) 

M2=M( 1.2) 

CORR=Wl 1.2)*FUN(M2) 

GO    TO    15 
11    M1=M(  II,  1) 

M2=M( 11,2) 

C0RR=WCI1,1)*FUN(M1)  +  W(  1 1 , 2 )*FUN { M2 ) 
15    M3=M2-1 

CO    21    K0=M1.M3 

011=FUN(K0)+FUN(KC+1) 

B=8+011 
21     CONTINUE 

R0=QU1)*(DEL*:B+C0RR) 

ROQ=ROO+RO 

E=.0 
7    CONTINUE 

ROOO=RC0*PP 

R000=(1.0/RATT )»RQOC 
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WRITEt6«40)RHO,ROOC,RQQO 
40  F0RMATC3D25.9) 

WRITEt7,204)RQQ0,PQCQ 

204  FORMATC 2F8.5) 
RQQ=0.0 

S  =  .0 

SY=-5.0 

IFtMM.EO.20)  GO  TO  67 

RHC=RHO+.05 

GO  TO  30 
67  RF0=RHO+.03 
30  CONTINUE   , 

205  CONTINUE 
STOP 
END 


C  TRANSFER    CHARACTERISTICS    OF    ODD    LEVEL    QUANTIZER 

C  FOLLOWED    BY    A    SQUARE    LAW    DETECTOR 

C  N=NUMBER    OF    STEPS 

C  VAR=NORMALIZED     STANDARD    DEVIATION    OF    X 

C  RHC=CORRELATICN    COEFICIENT    OF    X 

C  DEL=INTEGRATIOW    STEP 

C  XF=X(N) 

C  XC  =  X ( 0 ) 

IMPLICIT    REAL*8(  A-H,0-Z) 

DIMENSION    SYS(5001),Q(34),X(34),FUN(5  001)tW(18,2) 
lfM(18.2  ) 
C  PROBLEM    PARAMETERS 

N=17 

CO    205    NN=1,2 

READt5.200) VAR 
200  F0RMATIF4.1) 

WRITEt6.888)VAR 
888  FORMAT( '  • ,F4. 1) 
C      ALGORITHM  PARAMETERS 

RHC=0.0 

DEL=.001 

RQQ=0.0 

XF=1000. 

X0=-1000. 

SY=-5.0 

SS=0.0 

S=0.0 

B=0.0 

STP=(2./(N-1.) ) /VAR 

PP=. 3989423  . 

Nl=(N-l)/2 

L=DABS( SY) /DEL 
C      QUANTIZER   CHARACTERISTICS 

DO  81  1=1, N 

Ot  I)=( I-lN+1.)  /2.0)*STP 

Xf  I)  =  ( I-N/2.0)*STP 
81  CONTINUE 
C      LIMITS  OF  INTEGRATION 

DO  97  J  =  1,L 

SYS(J) =SY 

SY=SY+DEL 
97    CONTINUE 

DO    91    IR=1,N1 

IFtIR.ME.1)    GO    TO    92 

M(l,l )=1 

K3=tX(  D  +  5.0)  /DEL+1. 

W0=X(1)-SYS(K3) 

Mtl*2)=K3 

W( 1,2)=W0 

GO   TO    91 
92    PtIR.l)=VlTR-l ,2)+l 

W(  IR,1)=DEL-W0 

K4=lXt IR)+5.0)/DEL+l. 
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99 

101 

100 

98 
102 


W0=X(1P )-SYS(K4> 

M(  IR,2)=K4 

W(  IR,2  )=WO 
91    CONTINUE 

FINAL    POINT 

DO   98     1=1, N 

IF( I.NF.ll    GO    TO    99 

E2=X0/1.414 

E1=XU J/1.414 

GO    TO    100 

IFf.I.NE.N)    GO    TO    101 

E2  =  E1 

E1=XF/1.414 

GO    TO    100 

E2  =  E1 

F1=X{ I) /1.414 

R2=QERF<  El  )-DERF(E2) 

RR=0( I )**4*R2 

SS=SS+RR 

CONTINUE 

RATT=SS/2. 0 

RAT=DSQRT(RATT ) 

WRITE(6,102)RATT,RAT 

F0RMAT12F15.9) 

SY=-5.0 

DO    30    MM =1,21 

FORM    FUNCTION    LAMNA 

A=2.0*(l  -0-RH0**2) 

AA=DSQRT(A) 

DO    2    J=1,L 

DC    1    1=1, N 

IF(I.NE.l)    GO    TO    5 

E2=(X0-RH0*SY) /AA 

E1=(X(I )-RHO*SY)/AA 

GO    TO    25 
5     IF(I.NE.N)    GO    TO    3 

E2  =  E1 

El  =  (  XF-RHO^SY)  /AA 

GO    TO    25 
3    E2=E1 

E1=(X( I }-RHO*SY)/AA 
25    R1=D5RF(E1)-DERF(E2) 

R=0( I )**2*R1 

S=S+R 

1  CCNTINUE 

FUNC J)=( S/2.0J*DEXP< (SY**2>/(-2.0) ) 

S=.0 

SY=SY+DEL 

2  CCNTINUE 

MAIN    INTEGRATION 

CO    7    11=1, Nl 
44    IF( Il.NE.l)    GO    TO    11 

M1=M{ 1, 1) 

K2=M(1,2J 

CORR=W( 1,2)*FUN(M2) 

GO    TO    15 
11     M1=M(I1,1) 

M2=M( 11,2) 

CORR=W( I1,1)*FUN(M1)+W(I1,2)*FUN(M2) 
15    N3=M2-1 

DO    21    K0=M1,M3 

011=FUN(K0)+FUN(K0+1) 

B=B+Q11 
21  CONTINUE 

R0=0( II )**2*(DEL*B+C0RR) 

P0Q=FOQ+R0 

B=.0 
7    CONTINUE 

RQOO=RQQ-~-PP 

R000={ l.C/RATTJ*ROOO 

WRITE(6,40)RH0,RCQC,P000 
40    F0RMAT13D25.9) 
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WRITE(7,230)ROOO 

230 

FGRMAT(F8.5) 

ROQ=0.0 

S=.0 

SY=-5.0 

IFMM.EQ.20)  GO 

RHC=RH0+.05 

GO  TO  30 

67 

RH0=RHO+.03 

30 

CONTINUE 

205 

CONTINUE 

STOP 

END 

TO  67 


C      TRANSFFR  CHARACTERISTICS  OF  EVEN  LEVEL  QUANTIZER 

C      FOLLOWED  BY  A  SQUARE  LAW  DETECTOR 

C      N=NUMBER  OF  STEPS 

C      VAR=NORMALIZED  STANDARD  DEVIATION  OF  X 

C      RH0=CORRELATICN  COEFICIENT  OF  X 

C      DEL=INTEGRATION  STEP 

C      XF=X(N) 

C      X0=X(0J 

INPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  SYS( 5001)  ,0(34) ,X(34)  ,FUN ( 500 1 )  , W ( 18 ,2 ) 
1.MC18.2) 
C      PROBLEM  PARAMETERS 

N=4 

DO  205  NN=1,3 

READ(5,20Q)VAR 
200  F0RMAT(F4.1) 

WRITE(6,383)VAR 
888    FORMAT! »     • ,F4. 1) 
C  ALGORITHM    PARAMETERS 

RHO=0.0 

DEL=.001 

ROQ=0.0 

XF=1000. 

X0=-1000. 

SY=-5.0 

SS  =  0.0 

S=0.0 

6=0.0 

STP=(2./(N-1.) )/VAR 

PP=. 3989423 

Nl=N/2 

L=CA8S(SY)/DEL+1 
C  QUANTIZER      CHARACTERISTICS 

DO    81    1  =  1,  N 
..  0(I)=(I-(N+1.)  /2.0)*STP 

X(I)  =  (I-N/2.0)*STP 
81    CONTINUE 
C  LIMITS    OF    INTEGRATION 

DO    97    J  =  1,L 

SYS(J)=SY 

SY=SY+DEL 
97    CONTINUE 

DO   91     IR=1,N1 

IFCIR.NE.l )  GO  TO  92 

M(  1,11  =  1 

K3  =  (XU )+5 .01/DEL+l. 

WO=X( 1)-SYS(K3) 

Ml  1,2)  =  K3 

W(1,2)=WC 

GO  TO  91 
92  MUR.1)  =  M(  IR-1,21  +  1 

W(IR,1)=DEL-W0 

K4=(X(  IR1+5.0)  /DEL+1. 

W0  =  X(IR)-SYS(!<4) 

M(IR,2)=K4 
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W(  IR,2)=W0 

91  CONTINUE 

DO    91    IR=1,N1 

IF(IR.NE.l)    GO    TO    92 

M(l,l )=1 

DO    93    K3=l  ,L 

WO=SYS(K3)-X(l) 

IF(WO.GT.O.O)    GO    TO    94 

93  CONTINUE 

94  M(1,2)=K3-1 
W(l,2) =DEL-WO 
GO    TO    91 

92  M( IR,1)=M( IR-1,2)+1 
W(IR,l)=WO 

DO    95    K4=1,L 
W0=SYS(K4)-X(IR) 
IF(WO.GT.O.O)     GO    TO    96 

95  CONTINUE 

96  M(IR,2)=K4-1 
W(  IR  ,2)=DEL-W0 

91    CONTINUE 

FINAL    POINT 

DO    9  8    I  =1,N 

IF(I.NE.l)     GO    TO    99 

E2=X0/1.414 

E1  =  XU J/1.414 

GO    TO    100 
99    IF(I.NE.N)     GO    TO    101 

E2  =  E1 

E1=XP/1.414 

GO    TO    1J0 

101  E2=E1 

E1=X( I J/1.414 
100    R2=DERF(E1)-DERF(E2) 
RR=Q(  I  )*-*4*R2 
SS=SS+RR 
98    CONTINUE 

RATT=SS/2.0 

RAT=DSQRT(RATT) 

WRITE(6.102)RATT,RAT 

102  F0RMAT(2F15.9) 
SY=-5.0 

DO    30    KM =1,21 
FORM    FUNCTION    LAMNA 
A=2.0*(1.0-RHO**2) 
AA=DSORT(A) 
CO    2    J=1,L 
DO    1    1=1, N 
IF(l.NE.l)     GO    TO    5 
E2=(X0-RH0^SY)/AA 
E1=(X( I )-RHO*SY) /AA 
GO    TO    2  5 
5     IF (I .NE.K)     GO    TO    3 
E2=E1 

E1=(XF-RH0*SY)/AA 
GC    TO    25 
3    E2=E1 

E1=(X( I )-RHO*SY)/AA 
25    R1=DERF(E1)-DERF(E2  ) 
R=Q(I )**2*R1 
S  =  S+R 

1  CONTINUE 

FUN(J)  =  (  S/2.0)*DEXP(lSY**2)/(-2.0)  ) 

S=.0 

SY=SY+DEL 

2  CONTINUE 

^AIN    INTEGRATION 
DO    7    I  1  =  1, Nl 
44     IF(U.NE.l)    GO    TO    11 
M1=M(1 ,1  ) 
M2=M( 1,2) 
CORR=W( 1,21*FUN(M2) 
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GO  TO  15 
11  M1=M( 11,1) 

M2=M(I1,2) 

CORR=Wl II ,1)* FUN {Ml ) +W ( I 1 , 2 )*FUN { M2 ) 
15  M3=M2-1 

DO  21  K0=M1,M3 

011=FUN(KO)+FUN(KO+l) 

B=B+011 
21     CONTINUE 

R0=0U1)**2*(DEL*B  +  C0RR) 

ROO=RQO+RQ 

B=.0 
7  CONTINUE 

RQQO=ROQ*PP 

RGOC=(1.0/RATT )*RQQQ 

WRITE(6,40)RHO,RCCC,RQQO 
40  F0RMAT(3D25.9) 

WRITE(7,230)R0QC 
230  FORMAT(F8.5) 

ROQ=0.0 

S  =  .0 

SY=-5.0 

IF(MM.E0.20)  GO  TO  67 

RHO=RHO+.05 

GO  TO  30 
67  RHC=RH0+.03 
30  CONTINUE 
205  CONTINUE 

STOP 

END 


C  EXPECTED    VALUE    OF 

C  N=NUMBER    OF    STEPS 

C  XF=X{N) 

C  XO=X{0) 

IMPLICIT 

DIMENSION 

XF=100.0 

XC=-100. 

READ{5,71)N 
71     F0RMATU2) 

S=.0- 

VARR=.02 

CO    8    L=l ,440 

VAR=DSORT{ VARR) 

STP={2./(N-1.) )/VAR 

CO    1     1=1, N 

XU)  =  (  I-N/2.  )*STP 

0{I)=( I-(N+1.)/2.)*STP 
•  IFCI.NE.l)     GO    TO    5 

E2=X0/1.414 

E1=X(I)/1.414 

GO    TO    6 

5  IF(I.NE.N)     GO    TO    3 
E2  =  E1 
E1=XF/1.414 

GO    TO    6 
3    E2=E1 

E1  =  XU)  /1.414 

6  R1=DERF( E1)-DERF(E2) 
R=Q( I )**2*R1 

S  =  S+R 
1    CONTINUE 
QARR=VARR*S 
Y(L)=OARR 
V(L)=VARR 
SIG(L)=S/2.0 
VARR=VARR+.01 
S=0.0 


OUTPUT    AND    DERIVATIVE       I.F.     SAMPLING 


REAL*8[£-H,0-Z) 
0(34) ,X(34) , Y(440) , SI G 1440) ,V{440) ,Z(440) 
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8    CONTINUE 

CALL    DDET5( . 01 , Y,Z ,440, IER ) 

WRITE(6,82)(V(L),Y(L),Z(L),SIG(L) ,L= 1,440) 
82    FCRMATC     •  ,4F15.9) 
STOP 
END 


DEGRADATION    FACTOR    CF    IDEAL    FILTER 

IMPLICIT    REAL"8(  A-H,C-Z) 

REAL*8     ITITLE(12)/«  SWETT1,!!*'  •/ 

REAL*4    X3,F5,F6,F7,F8,F9,F10 

REAL*8    TITLE( 10J/10*'     •/ 

REAL*4     LABEL/4H  / 

DIMENSION    Z(43)  ,F2(43)  ,  W  I  (43  )  ,  Y  (43")  ,  DELY  (  43  )  ,  B  (21  ) , 
1F(9,31),F5(31),F6(31),F7(31),F8(31),F9(31),F10(31), 
2AP1  (31),  AP2(  31  ),AP3(31),AP  4(31),  AP  5(31),  AP  6(31), 
3F  12(31)  ,F13(31)  ,SB(21)  ,X3(31) 

DO  98  1=1,22 

98  READ(5,99) WI (I ) 

99  F0RMAT(F7.1) 
DO  55  MM=1,6 
READ(5,999)VAR1,VAR2,VAR3,VAR4,VAR5 

999  FORMAT ( 5 F4.1) 
CO  5  M  =  2,6 
DO  1  1=1,22 
READ(5,2)F2(I ) 

2  F0RMATIF8.5) 
1    CONTINUE 

READ(5,47)AA,E,C,D 

47  F0RMAT(4F15.7) 
A=AA-C**2 
WRITE(6,48)A,E,C,D 

48  FORMAT( •     • ,4F15.7) 
Z(21)=.98 

Z(22)  =  l  .0 

S=0.0 

DO    3    1=1,20 

Z(I)=S 

S=S+.05 

3  CONTINUE 

;  FIT    POLYNOMIAL 

CALL    LS0PL2(22,-16,Z,F2,WI ,Y, DELY, B, SB, TITLE) 

E(1)=0.0 

bT=100. 0 

XX=.01 

DO    35    1=1,31- 

N=WT/XX-1.0 

RUN=0.0 
.  DO    4    J=l ,N 

IF(J.GT.650)    GO    TO    49 

T=3.1416*J*XX 

X7=DSIN(T)/T 

X=DABS( X7) 

IF(X.LT.O.l)    GO    TO    3332 

PUN=B(  1)+B(2)*X+3(3)*X**2  +  B(4)*X#*3  +  B(5)*X**4- 
1+B(6)*X**5+B(7)*X(6  )  +  B( 8 )*X**7+B { 9)*X< 8  ) +B  (  1  0  )*X**9 
2+B(ll )*X**10+3(12)*X**ll+B{13)*X**12+3< 14)*X**13 
3+B( 15)*X**14+3{16)#X**15+B(17)*X**16 

GO    TO    3333 

3332  PUN=(F2(3)-F2(1 ) )*X*19.0 

3333  RUN=RUN+PUN 

4  CONTINUE 

49  A1=A*XX 
A2=2.0*XX*E*RUN 
F4=2.0*(  Al+£2)/D**2 
F(M,I)=DSQRT(F4) 
IF(I.NE.l)    GO    TO    66 
XX=XX+-09 

GO    TO    35 
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66    XX=XX+.l 
35    CCNTINUE 
5    CONTINUE 
XX=.l 

F(l,l)=1.0 
X3(l)=0.0 
I  NO    QUANTIZER 

DO    7    1=2,31 
N=WT/XX-1.0 
R  =  0.0 
DO    8    J  = 
T=3.141 
S=(DSIN 
R=R+S 
8    CONTINU 
YY=XX+2 
F(1,I)= 
X3(I)=X 
XX=XX+. 
7    CONTINU 
WRITE<6 
9  87    FORMAT ( 
WRITE(6 
1F(6,I ) , 
31     FORMAT { 
WRITE16 
888    FORMATt 
X=0.0 
READ(5, 

88  FORMAT{ 
WRITEC6 

89  FORMAT( 
;      APPROXI 

DO    297 

RR=F(1, 

SS=X*( 2 

SSS=X*( 

APP=RR+ 

APK  J)  = 

APP=RR+ 

AP2( J)= 

APP=RR+ 

AP3(J)= 

APP=RR+ 

AP4(J)= 

APP=RR+ 

AP5(J)= 

X=X+.l 
297    CONTINU 

WRITE(6 

1,AP5( J) 

1100    FORMAT ( 

CO    18    K 

F5(K)=F 

F6(K)=F 

F7(K)=F 

F8(K)=F 

F9(K)=F 
18    F10(K)= 

CALL    DR 
16,1, LAS 

CALL    DR 
16, 1,LAS 

CALL    DR 
16,1, LAS 

CALL    DR 
16,1, LAS 

CALL    DR 
16, 1,LAS 

CALL    DR 

1,6,1, LA 

55    CONTINU 


CASE 


1,N 

6*J*XX 
(Tl/T)*  = 


.0*XX*R 

DSORT(YY) 

X 

1 

r 

t?87) 

•  1«  ) 

,31)  (X3(  I 

1=1,31) 

'     «,7F12. 

,888) 

88)N 

12) 

,89)N 

'     ' ,12) 

MATION 

J=l,31 

J)**2 

.0/CN-1J ) 

1 .0/360  .G 

SS/VAR1** 

DSQRTIAPP 

SS/VAR2** 

DSORT( APP 

SS/VAR3** 

DSORTIAPP 

SS/VAR4*- 

DSORT( APP 

SS/VAR5+* 

DSORTIAPP 


),F(1,I),F(2,I),F{3,I),F(4,I),F(5,I), 

7} 

) 


, 1100)( X3 

, J  =  l,31  ) 

'     «,7F12. 

=  1,31 

(l.K) 

(2.K) 

(3,K) 

(4,K) 

(5,K) 

Fi6,K) 

AM31.X3. 

T) 

AW(31,X3, 

T) 

AW(31,X3, 

T) 

AW(31,X3, 

T) 

AW131  ,X3, 

T) 

AW(31  ,X3, 

ST) 

E 


**2/6.0 
)*«2.0/{N-lJ  )**4 

2+SSS/VARl**4 

) 

2+SSS/VAR2**4 

) 

2+SSS/VAR3**4 

) 

2+SS5/VAR4**4 

) 

2+SSS/VAR5**4 

) 


(J),F( 1, JJ ,AP1( J) ,AP2(J) ,AP3( J) ,AP4(J) 
7) 


F5, 1, 0, LABEL, IT IT LE, 0.5, 1.0, 0,0, 2, 2, 6, 
F6, 2,0, LAB EL, IT ITLE,0 .5 , 1.0,0,0, 2, 2, 6, 

F7,2,0,LABEL,ITITLE,0.5,1. CO, 0,2, 2,  6, 
F 8,  2,0,  LABEL, I  TITLE, 0.5, 1.0, 0,0,2,2,6, 
F 9, 2,0, LAB EL, I TITLE, J. 5, 1.0, 0,0, 2, 2, 6, 
F10,3,0,LABL5, ITITLE, 0.5, 1.0,0,0,2, 2, 6 
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STOP 
END 


DEGRADATION  FACTOR  OF  SINGLE  TUNED  FILTER 

IMPLICIT  REAL*3( A-H.O-Z) 

REAL*4  X3,F5,F6,F7,F8,F9,F10 

REAL*8  ITITLE112)/'    SWETT',11*'  '/ 

REAL*8  TITLEtl J)/10*»  *J 

REAL*4  LABEL/4H     / 

DIMENSION  Z(43 ) , F2 ( 43 ) , W I { 43 ) , Y ( 43)  ,DELY( 43)*B(21) , 
1F(9,31),F5(31),F6(31),F7(31),F8(31),F9{31),F10(31), 
2AP1(31),AP2(31),AP3(31),AP4(31),AP5{31),AP6(31), 
3F12C31)  ,F13(31 ),SB(21) ,X3(31) 

DO  98  1=1,22 

98  READ{5.99)WI(I  ) 

99  FCRMAT(F7.1) 
DO  55  MM=1,6 
REA0(5,999)VAR1,VAR2,VAR3,VAR4,VAR5 

999    F0RMAT15F4.1) 
DO    5    M  =  2,6 
CO    1    1=1,22 
READ(5,2)F2(I ) 

2  FORMAT (F8. 5) 
1     CONTINUE 

READ( 5,47)AA,E,C,D 

47  FORMAT (4F 15. 7) 
A=AA-C"*2 
WRITE(6,48) A,E,C,D 

48  FORMAT!'     S4F15.7) 
Z(21)=.98 
Z(22)=1.0 

S=0.0 

DO    3    1=1,20 

Z(  I)  =  S 

S=S+.05 

3  CONTINUE 

FIT    POLYNOMIAL 

CALL    LSCPL2(22,-16,Zt F2 , W I , Y ,DELY, B, SB, TITL E ) 

B(1)=0.0 

WT=10.0 

XX=.01 

DO    35    I =1,31 

N=WT/XX-1.0 

RUN=0.0 

DO    4    J=1,N 

IF(J.GE.2CO).   GO    TO    4 

T=3.1416*XX=i'J 

X=DEXP(-T) 

IF(X.LT.O.l)    GO    TO    3332 

PUN  =  B(1)+B(2);f-X+E(3  )*X**2+3 (4 )*X**3+B ( 5  )*X**4 
l  +  8(6)*X**5  +  B(7)-X**6+Bl8)*X*~7+B(9)*X*:fa8  +  B(10)*X**9 
2+B(ll )*X**10+B(12)*X**11+B( 13)*X**12+B( 14)*X**13 
3+B ( 1 5 ) *  X** 14+B (16)*  X** 15  +  B (17) *X**16 

GO    TO    3333 

3332  PUN=(F2 (3)-F2( 1 ) )*X*10.0 

3333  RUN=RUN+PUN 

4  CONTINUE 
Al=A=r-XX*3.1416 
A2=2.0*3.1416*XX*E*RUN 
F4=2.0-*(A1  +  A2) /D**2 
F(M,I)=DSQRT(F4) 
IFCI.NE.l)  GO  TO  66 
XX=XX+.09 

GC  TO  35 
66  XX=XX+.l 
35  CONTINUE 

5  CONTINUE 
XX  =  .l 
F(l,l)=1.0 
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X3(l)=0.0 

NO  QUANTIZER  CASE 

DO  7  1=2,31 

R=3.1416*XX 

YY=R/DTANH(R) 

F(1,I)=DSQRT(YY) 

X3(I)=XX 

XX=XX+.l 
7  CCNTINUE 

WRITE(6,<587) 
987  FORMAT (• 1*  ) 

WRITF(6,31)(X3(  I)  ,FU,I)  ,F(2,I)  ,F(3,I),F(4,  I),F{5,  I), 
1F(6,I ) , 1=1,31) 
31  FORMAT  I ■  »,7F12.7) 

WRITE(6,888) 
888  FORMATt*  •,/////) 

X=0.0 

RFAD(5,88)N 

88  FORMAT(I2) 
WRITE(6,89)N 

89  FORMAT( «  • ,12) 
;      APPROXIMATION 

DO  297  J=l,31 

RR=F(1,J)**2 

SS=X*(2.0/(N-1  )  )**2/6.0 

SSS=X*{1.0/360.0)*(2.0/(N-1)  )**4 

APP=RR+SS/VAR 1**2+ SSS/VAR 1**4 

APl(J)=DSORT(APP) 

APP=RR+SS/VAR2**2+SSS/VAR2*-:*4 

AP2( J)=DSORT(APP) 

APP=RR+SS/VAR3*:*2+SSS/VAR3**4 

AP3( J)=DSQRT(APP) 

APP=RR+SS/VAR4**2+SSS/VAR4**4 

AP4(J) =DSQRT(APPl 

APP=RR+SS/VAR5**2+SSS/VAR5**4 

AP5(J)=DSQRT(APP) 

X=X+.l 
297  CONTINUE 

WRITE{6,1100HX3(J),F(1,J),AP1(J),AP2(J),AP3(J),AP4(J) 
1,AP5( J)  ,J=1,31) 
1100  FORMAT ('  «,7F12.7J 

DO  18  K=l,31 

F5(K)  =  F{  1,K) 

F6(K  )=F(2,K) 

F7(K)=P(3,K) 

F8(K)  =  F(4,K) 

F9(K)=F(5,K) 
18  F10(K)=FI6,K) 

CALL  DP.  A  W(  3  1,X3,F5,  1,0,  LABEL,  I  TITLE,  0.5, 1.0,  0,0, 2,2,  8, 
16.1, LAST ) 

CALL  DRAWC31  ,X3 , F6 , 2 ,0 , LABEL , IT  IT LE ,0 .5 , 1 . 0, 0 , 0 ,2 , 2 ,6 , 
16,1, LAST) 

"  CALL  DRAW131, X3,F7, 2,0, LABEL, ITITLE, 0.5, 1.0,0,0,2,2,6, 
16,1, LAST) 

CALL  DRAW(31,X3,F8,2, C, LABEL, ITITLE ,0.5,1.0,0,0,2,2,6, 
16,1, LAST) 

CALL  DRA  W(  31  f  X3 1  F9*  2t  0  f  LABEL  tl  TITLE  1 0.5  ,1.0,0,0,2,2,6, 
16,1, LAST ) 

CALL  DRAW(31,X3,F10,3,3,LABL5,ITITLE,0.5,1.0,0,0,2,2t6 
1,6.1, LAST) 
55  CCNTINUE 

STOP 

END 


DEGRADATICN  FACTOR  OF  GAUSSIAN  FILTER 

IMPLICIT  REAL*8( A-H.C-Z1 

REAL*4  X3,F5,F6,F7,F8,F9,F10,F11 

REAL*8  TITLEdO)/^*1  '/ 

RFAL*4  LABEL/4H     / 


lHO 


58 
99 


999 


2 

1 

47 


48 


3332 

3333 

4 


66 

35 
5 


8 


DIMENSION    Z(43 
1F(9,31  )  ,F5  (31) 
2APK31)  ,AP2(31 
3F12I31), F13131 
DO    98    1=1,22 
READ(5,99)WI (I 
F0RMAKF7.1  ) 
DO    55    MM=1,2 
READ(  5,999)VAR 
F0RMAT13F4.1  ) 
DO    5    M  =  2,4 
DO    1    1=1,22 
READ(5,2)F2( I  } 
FORMATIF8. 5) 
CONTINUE 
READ(5,47) AA,E 
F0RMAT(4F15.7) 
A=AA-C**2 
WRITE(6,48)A,E 
FORMATC     «,4F1 
Z(21)=.98 
Z(22)=1.0 
S=0.0 

DO   3    1=1,20 
Z(I)=S 
S=S+.05 
CONTINUE 
FIT    POLYNOMIAL 
CALL    LSQPL2(22 
U=2.773 
B(  1)  =  0.0 
WT=10.0 
XX=.Q1 

DO    35    1  =  1,30 
N=WT/XX-1.0 
RUN=3.0 
DO    4    J=1,N 
IF(J.GE.IOO)    G 
T=( (3.1416-XX* 
X=DEXP(-T) 
IF(X.LT.O.l)    G 
PUN=B( 1J+B12)* 
1+B(6)*X**5+B(7 
2+B(ll)*X**10+3 
3+B(  15)*X**14+B 
GO   TO    3333 
PUN=(F2(3)-F2( 
RUN=RUN+PUN 
CONT  INUE 
A1=A*XX*1.505 
A2  =  2.0'-l:1.50  5irX 
F4=2.0*(A1  +  A2) 
F(M,I ) =DSORT(F 
IF(I.NF.l)    GO 
XX=XX+.09 
GO    TO    35 
XX=XX+. 1 
CONTINUE 
CONTINUE 
XX=.l 

F  (  1 , 1 )  =  1  .  0 
NO    OUANTIRZER 
CO    7    1=2,30 
N  =  WT/XX-1.  3 
R  =  0.0 
X3(1J=0.0 
S=0.0 

DO    8    J=1,N 
T=(2.3* (3.1416 
S  =  DEXP(-T) 
R  =  R+S 
CONTINUE 


), F2143) ,WI(43) ,Y(43) ,DELY(43) ,B(21) 
,F6(31),F7(31),F8(31),F9(31),F10(31) 

) 
) 


AP3(31) ,AP4(31) , APS (31) , AP6(31}, 
SB(21) ,X3(31) ,F11( 31) 


1, VAR2.VAR3 


t  CD 


.CD 
5.7) 


,-16,Z,F2,WI,Y,DELY,B,SB,TITLE) 


0    TO    4 
J)**2)/U 

0  TO    2332 

X+  B ( 3 ) *  X** 2+  B 1 4 ) *  X**3  +  B ( 5 ) *X**4 
)*X**6+B(  8)=*X**7+B(9)*X**8  +  B(  10)*X**9 
(12)*X**11+B113)*X**12+B(  14)=pX~*13 
(  16)*X-**15+B(17)=*X**16 

1  )  )*X*10.0 


X~E*RUN 

/C**2 

4) 

TO    66 


CASE 


*XX*J)**2)/U 
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987 

31 

888 


88 


89 


V=1.505* 
FC1.  I)  =  D 
X3(I)=XX 
XX=XX+.l 

CONT  INUE 

WRITE (6, 

FORMAT! * 

WRITE(6, 

FORMAT! • 

WRITE16, 

FORMAT  {• 

X=0.0 

READ!5,8 

FORMAT! I 

WRITE(6, 

FORMAT! • 

APPROXIM 

DO    297    J 

RR=F!1, J 

SS=X*(2. 

SSS=X*(  1 

APP=RR+S 

AP1!J)=D 

APP=RR+S 

AP2! J)=D 

APP=RR+S 

AP3!J)=D 

X=X+.l 

CONTINUE 

WRITE16, 

DO    18    K= 

F5!K)=F! 

F6!K)=F! 

F7!K)=F( 

F8!K)  =  F( 

F9(K)=AP 

F10!K)=A 

F11(K)=A 

CALL    DRA 

16,1,  LAST 
CALL    DRA 

16,1, LAST 
CALL    DRA 

16,1, LAST 
CALL    DRA 

16, 1,LAST 
CALL    DRA 

16,1, LAST 
CALL    DRA 

1,6,1, LAS 
CALL    DRA 

1  ,6,1, LAS 
55    CONTINUE 
STOP 
END 


297 


19 
20 
21 


18 


XX+3.01*XX*I 
SORTIV  ) 


987) 
1"  ) 
31)(X3II),F(1,I),F(2,I),F!3,I),F{4,I},I=1,3  0) 

• ,5F12.7) 
888) 

'  ,/////) 

8)N 
2) 
89  )N 

'  ,12) 
ATICN 
=  1,30 
)*-*2 

O/CN-1 ) )**2 
.0/360.0)*! 
S/VAR1**2+S 
SORT(APP) 
S/VAR2**2+S 
SORT(APP) 
S/VAR3**2+S 
SQRT(APP) 


/6.0 

2.0/(N-l)  )**4 
SS/VAR1**4 

SS/VAR2**4 

SS/VAR3**4 


31)  (X3!J),F 

1,30 
1,K) 

2,K) 

3,K) 

4,K) 

1!K) 

P2(K) 

P3!K) 

W(30,X3,  F5, 

)) 

W!30,X3, F6, 

) 

lrJ(30,X3,F7, 

3 

W(30,X3,F8, 

) 

W!30,X3,F9, 

) 

W(30,X3,F10 

T) 

W! 30,X3,F11 

T) 


!1,J),AP1!J),AP2!J},AP3(J),J=1,30) 


1,0, LAB  EL, IT ITLE, 0.5, 0.5, 0,0,  2, 2,  6, 
2,0, LAB  EL, I  TITLE, 0. 5, 1. 0,0, 0,2, 2 1 6, 
2,0,  LAB EL, I  TITLE, 0.5 ,1.0,0,0,2,2,6, 
2,0, LABEL, IT ITLE, 0.5, 1.0, 0,0, 2, 2, 6, 
2,  3,  LAB  EL,  I  TITLE,  0.5,1.  CO,  0,2,  2, 6, 
,2,4,LABEL,ITITLE,0.5,1.C,0,0,2,2,6 
,3, 5,  LABEL, I  TITLE, 0.5, 1.0, 0,0, 2, 2, 6 


DEGRADATICN    FACTOR    C=    BUTTERWORTH 
IMPLICIT    REAL*8( A-H,0-Z) 
REAL*4    X3,  F5,F6,  F7,F3,F9,F10 
REAL*8    ITITLE112)/' 
REAL*8    TITLE! 10)/10*' 
REAL*4    L£BEL/4H  / 


FILTER 


SWETT* 
«/ 


11* 


DIMENSION  Z(43) 
1F(9,31) ,F5(31) , 
2APK31)  ,AP2!31  ) 
3F12!31)  ,F13(31i 

CO    98     1=1,22 

98  READ(5,99)WI  (I  ) 

99  F0RMAT!F7.1) 


F2(43) 
F6(  31) 


,WI(43) 
F7!31) 


,Y!43) 

F8!31) 


,DELY(43),B!21 ) 
,F9!31)  ,F10(31J 


AP3(31),AP4(31),AP5(3i),AP6(31) 
, SB(21)  ,X3(31) 
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CO    55    MM=1,6 

RE AD (5. 999) VAR1 , VAR2  ,VAR3 » VAR4 t VAR5 
999    FORMAT ( 5F4.1) 
CO    5    M=2,6 
DO    1    1=1,22 
READ(5,2)F2(  I  ) 

2  FOPMAT(F8.5) 
1    CONTINUE 

RFAD(5,47)AA,E,C,D 

47  F0PMAT(4F15.7) 
A=AA-C:':*2 
WRITE(6,48,)A,E,C,D 

48  FORMAT! •     • ,4F15.7) 
Z(21)=.98 
Z(22)-1.0 

S  =  0.0 

DO    3    1=1,20 

Z(I)=S 

S=S+.35 

3  CONTINUE 

FIT    POLYNOMIAL 

CALL    LSQPL2(22,-16, Z,F2,WI ,Y,DELY, B, SB, TITLE) 

B(1)=0.0 

WT=10.0 

XX=.01 

DO    35    1  =  1,31 

N=WT/XX-1.0 

RUN=0.0 

DO    4    J=1,N 

1F(J.GT.30J)    GO    TO    49 

T=1.414*DEXP(-J*3. 14  16*XX/1.414) 

X7=T*DC0S( J*3.14  16*XX/ 1.414-3.1416/4.0) 

X=DABS( X7) 

IF(X.LT.O.l)    GO    TO    3332 

PUN=B(1 )+B(2)*X+E(3 )*X**2+B ( 4 )*X**3+B( 5 )*X**4 
1+B{6)*X**5+3(7)*X**6+R<8)*X**7+B(9)*X**8+B(10)*X**9 
2  +  B ( 1 U  *X  ** 1 0  +B ( 1 2 ) *  X** 1 1  +  B  < 1 3  )  *  X**  1 2+B  (14)*  X**  1 3 
3+B(15)»X**14+BU6)*X#=n5+B(17}*X**16 

GO    TO    3333 

3332  PUN=(F213)-F2( 1) )*X*10.0 

3333  RUN=RUN+PUN 

4  CONTINUE 

49  A1=A*XX*1.48 

A2  =  2  .0*1  .48^- XX*  E*RUN 

F4=2.0*(A1+A2>/D**2 

F(M, I)=DSQRT(F4) 

IFd.NE.  1)     GO    TO    66 

XX=XX+.09 

GO    TO    35 

66    XX=XX+.  1 

35    CONTINUE 

5.  CCNTINUE 

XX=.l 

F(  1.1). =  1.0 

X3(l)=0.0 

U=1.414 

NO    QUANTIZER    CASE 

DO    7    1=2,31 

N=WT/XX-1.  0 

R  =  0.0 

DO    8    K=1,N 

0=3.1416*KiXX*U 

S=DEXP(-0)*(1.0+DSIN(0) ) 

R  =  R+S 
8    CONTINUE 

YY=1 .48TXX+2 .96*XX*R 

F( 1,I)=DS0RTIYY) 

X3( I )=XX 

XX=XX+.l 
7    CONTINUE 

WRITE(6,987) 
987    FCRMATC  1«  ) 
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WRITE (6, 31)  (X3(I)  ,F(1,I)  ,F(2,I)  ,F(3, I)>  F(4,  I), F(5,l), 
1F(6, I  ), 1  =  1,31) 
31    FGPMATC     «,7F12.7) 
WRITE(6,888) 
888    FORMAT!'     ',/////) 
X=0.0 
READ!5,88)N 

88  FORMATU2) 
WRITE!6,89)N 

89  FORMAT!  '     '  ,12) 
APPROXIMATION 
DO    297    J=1.31 
RR=F( 1,J)**2 
SS=X*!2.0/!N-1 )  )**2/6.0 

SSS=X* (1.0/360.0)*! 2. 0/(N-l) )**4 

APP=RR+SS/VAR1**2+SSS/VAR1#*4 

APK J)=QSORT(APP) 

APP=RR+SS/VAR2** 2+SSS/ VAR2**4 

AP2( J)=DSORT( APP) 

APP=RR+SS/VAR3**2+SSS/VAR3**4 

AP3(  J)-=DSQRTIAPP) 

APP=RR  +  SS/VAR4*:"2  +  SSS/VAR4**4 

AP4( J)=OSORT( APP) 

APP=RR+SS/VAR5**2+SSS/VAR5**4 

AP5(J)=DSQRT( APP> 

X=X+.l 
297  CONTINUE 

WRITE(6,1100)(X3(J),F(1,JV,AP1(J),AP2(J),AP3(J) ,AP4( J) 
lt-AP5(  J)  .0=1,31) 
1100  FORMAT! •  «,7F12.7) 

DC  18  K=l,31 

F5(K)=F(  1,K) 

F6(K)=F(2,K) 

F7(K)=F(3,K) 

F8(K)=F(4,K) 

F9(K)=F(5,K) 
18  F10(K)=F(6,K) 

CALL  DRAW! 3 1 , X3, F5 t 1 , 0 , LABE L , I TI TLE ,0.  5  , 1 .  C  ,  0  ,0  ,2  ,2  ,6  , 
16,1, LAST) 

CALL  DRAW! 31 , X3 , F6 , 2 , 0 , LABEL , IT  IT LE , 0 .5 ,1.0,0,0,2,2,6, 
16,1, LAST) 

CALL  DRAW!31,X3, F7, 2 , 0, LABEL , ITI TLE, 0.5, 1.0,0,0,2,2,6, 
16, 1,LAST) 

CALL  DR AW! 3 1,X 3, F 8, 2,  CUL ABEL, I  TITLE, 0.5  ,1.0,0,0,2,2,6, 
16,1, LAST) 

CALL  DRAW! 3 1 , X3 , F9 , 2, 0 , LABEL , I TI TLE , 0. 5 , 1 . 0  ,0  ,0 ,2 ,2,6, 
16.1, LAST  ) 

CALL  DRAW(31,X3,F10,3,0,LABL5,ITITLE,0.5,1.0,3,0,2,2,6 
1,6.1, LAST)  • 
55  CCNTINUE 

STOP 
.  END 


DEGRADATION    FACTOR.       RC    FILTERING    AND    SAMPLING 

IMPLICIT    REALMS! A-H,C-Z) 

REAL*8    ITITLEC 12) /•  SWETT«,11*«  •/ 

REAL^4    X3,F5,F6,F7, F3, F9,F10 

REAL*4    LABEL/4H  / 

REAL*8    TITLE!  10)/10*«     •/ 

DIMENSION    Z(43)  , F2 (43) , W  I  (43 ) , Y (43 ) , DELY ( 43  ) , 8( 2 1 ) , 
1F(9,31).F5(31)  ,F6!31),F7(31),F8(31),F9!31)  ,F10(31), 
2APK31),AP2(31),AP3f31),AP4(31),AP5(31J,AP6(31), 
3F12131).F13(31),SB(21),X3(31) 
4,FC(9,31),E(20) ,0(20) 

GAUSS  (E,VAR)=0. 3989423*  ( DEXP (- ( E/VAR )**2/2  .  )  ) /VAR 

DO    93    1=1,22 

98  READ(5,99)WI(I ) 

99  FCRMATIF7.1) 
DO    55    MM=1,6 
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144 


88 

89 

999 


2 
1 

47 


2000 


2001 


2500 
79 

48 


3332 
3333 


READ(5.8 
FGRMATU 
WRITEl 6, 
FORMAK  • 
READ15.9 
FORMAT! 5 
CO    5    M=2 
DO    1    1=1 
READ(5,2 
FORMAT (F 
CONTINUE 
READ15.4 
FORMAT (2 
STP  =  ( 2/1 
FF=0.0 
DER=0.0 
CALCULAT 
DO    79     1= 
0(I)=(I- 
F( I)=(I- 
I  F  (  I  .  N  E  . 
E3  =  EU ) 
E1=0.0 
E2  =  GAUSS 
GO    TO    2  5 
IFd.NE. 
E1=E2 
E2=0.0 
GO    TO    25 
E4  =  EU ) 
E1  =  E2 
E2=GAUSS 
PPEL=Q(I 
DER=DER+ 
WRITE(6, 
FORMAK1 
Z(21)=.9 
Z(22)=l. 
S=O.Q 
DO    3    1=1 
Z(I)=S 
S=S+.05 
CONTINUE 
FIT    POLY 
CALL'LSQ 
TRC=20.0 
XX=.01 
DO    35    1= 
NN=TRC/X 
RUN=0.0 
RRUN=0.0 
DO    4    J=l 
T=XX*J 
X=DEXP(- 
IFIX.LT. 
PUN=B(1 ) 
1+B(6)*X* 
2+9(ll)*X 
3+B(15)*X 
GO    TO    33 
PUN=F2(3 
RUN=RUN+ 
RRUN=RRU 
CONTINUE 
A1=A*XX/ 
A2=XX*A* 
A3=XX**2 
F41=A1+A 
F42=A1+A 
FCM.I J=D 
FC ( M ,  I  )  = 
IFd.NE. 


8)N 
2) 

89)N 
•  ,1 
99)  V 
F4.1 
,6 
.22 
)F2( 
8.5) 


2) 

AR1 , VAR2  .V AR3 , VAR4, VAR5 
) 


I  ) 


7JA.VAR 
F15.7) 
N-1.0)  ) 


ICN    OF    DERIVATIVE    OF    EXPECTED    VALUE 

l.N 

(N+1.0)/2.0)*STP 

N/2.0)*STP+FF 
1)     GO    TO    2000 


IE3.VAR) 

N)    GO    TO    2001 


00 


(E4.VAR) 
)*(E1-E2) 

PPEL 
48)A.VAR,DER 

•  .3F15.7) 
8 
0 

,20 


NOMIAL 

PL2(22 ,-16,Z,F2,WI ,Y,DELY,BtSBtTITLE) 


1.31 
X-1.0 


,NN 


T) 

0. 
+  B 
*5 


23 
)* 

PU 
N+ 


1)    GO    TO    3332 

(2)*X+E(3  )*X**2  +  B(4)*X*#3  +  B(5)*X=**4 

+B(7)*X**6+B(8)*X**7+B(9)*X**8+8(10)*X**9 

10+B(  12)~X-"*11+B(  13)^X-^*12  +  3(  14)*X**13 

14+3(16)*X**15+B(17)*X**16 

10.0*X 

N 
J*PUN 


2.0 

RUN 

*A*RRUN/TRC 

2-A3 

2 

S0RKF42 )/CER 

DSQRTl F41J /OER 

1)     GO    TO    66 
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XX=XX+.09 

GO    TO    35 

66    XX=XX+.l 

35    CONTINUE 

5    CONTINUE 

XX  =  .l 

F( 1,11=1.0 
YYY=1.0-(  1 
FCU,11=DS 
X3( 1)=3.0 
NO  QUANT  IZ 
DO  7  1=2,3 
YY=(XX/2.0 
Y4=(XX/2.0 
Y3=YY-(1. i 
F( 1,  I)=DSQ 
FCU.  I  )=ns 

X3(I) =XX 
XX=XX+.l 
7    CONTINUE 
WRITE  J  6,98 
987    FORMAT  (  •  1' 
WRlTE(6,31 

1F(6,I) ,1=1 
X=0.0 

APPROXIMAT 
DO  297  J=l 
RR  =  F(1,  J  )* 
SS=X*(2.0/ 
APP=RR+SS/ 
AP1(J}=DS0 
APP=RR+SS/ 
AP2(  JJ=DSO 
APP=RR+SS/ 
AP3( J) =DSO 
APP=RR+SS/ 
AP41J)=DSQ 
APP=RR+SS/ 
AP5( J)=DSO 
X=X+.l 
CONTINUE 
WRITE(6,88 
FORMAT! •  • 
WRITEC6, 11 

1  ,AP5( J) , J  = 
1100  FORMAT! '  • 
WRITE(6,  13 
FORMAT!'  • 
FORMAT! '  • 
DO  18  K=l, 
F5!K)=F!1, 
F6(K)=F(2, 
F7(K)=F(3, 
F8!K)=F!4, 
F9(K)  =  F(  5, 
18  F10(K)=F!6 
CALL    DRAW( 

16, 1,LAST) 
CALL    DRAW! 

16,1, LAST) 
CALL  DRAW! 

16,1, LAST) 
CALL  DRAW! 

16,1, LAST ) 
CALL  DRAWl 

16,1, LAST) 
CALL  DRAW! 

1,6,1, LAST) 
55  CONTINUE 
STOP 
END 


297 

888 


13 
31 


.O-DEXPt-TRC) )/TRC 
QRTIYYY) 

ER    CASE 

1 

) /DTANH(XX/2.0) 

l/DSINHJXX/2.0) 

-CEXP(-TRC)  )*Y4**2/TRC 

RT(YY) 

GRTIY3  ) 


7) 
) 

)  !X3!  I)  ,F!1,  I),F(2,  I),F!3,  I),F{4,  I),F(5,I), 
,31) 

ICN 

,31 

*2 

(N-1.0) )**2/24.0 

VAR1**2 
RT( APP  ) 
VAR2^*2 
RT( APP) 
VAR3*~2 
RT( AFP) 
VAR4-*2 
RT(APP) 
VAR5**2 
RT( APP) 


8) 

,/////) 

00MX3!J),C(1,J),AP1(J),AP2(J),AP3(J),AP4(J) 

1,31  ) 

,7F12.7) 

) 

,/////) 

,7F12.7) 

31 

K> 

K) 

K) 

K) 

K) 

,K) 

3 1,X3,F6, 1,0, LABEL, I  TITLE, 0.5, 0.0, 0,0, 1,2, 6, 

31,X3,F5,2,C,LABEL, ITITLE,0.5,C.C,0,0,1,2,6, 

3  1,X3, F 7, 2,0, LAB EL, I  TITLE, 0.5,0. 5,  0,0, 2, 2, 6, 

3 1,X3,F 8, 2,0, LAB EL, I  TITLE  ,0.5 ,0.5,0,0,2 ,2,6, 

3 1,X3,F9,  2,0,  LAB  EL,  I  TITLE,  0.5,  3.  5  ,-0,0,  2,  2,  6, 

31,X3,F10,3,0,LABEL,ITITLE,0.5,C.5,0,0,2,2,6 
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C  TRANSFER  C 
C  N=NUM6ER  0 
C  VAR=NORMAL 
C  RHO=CORREL 
C  DEL=INTEGR 
C  XF=X(N) 
C  XO=X(0) 
C      F=OFFSET 

IMPLICIT  R 
DIMENSION 
C      PROBLEM  PA 
N=13 

WRITEI6.66 
666  FORMAT( •  ' 
VAR=0.2 
C      ALGORITHM 
RH0=-0.999 
DEL=.001 
F=-0.5 
ROO=0.0 
XF=10.0D  1 
XO=-10.0D 
SY=-5.0 
SS=0.0 
S  =  0.0 
B=0.0 

STP=(2./  IN 
PP=. 398942 
L=10.0/DEL 
C  QUANTIZER 
DO  81  1=1, 
0(I)  =  (  I-CN 

xm=(  (i-n 

81  CONTINUE 
C      LIMITS  OF 
Nl  =  l 
N2=N 
N3=N-1 
DO  911  IRR 
IFIXl IRR). 
IFIXIIRR). 
GO  TO  911 
10CO  N1=IRR+1 
GO  TO  911 

911  CONTINUE 
GO  TO  912 

1001  N2=IRR 

912  DO  91  IR=N 
IF(  IR.NE.N 
MIIR,1)=1 
K3=IXI IR}  + 

..  WO=XIIR)-( 
Ml  IR,2)=K3 
W(  IR  ,2.)  =  W0 
GO  TO  91 
922  IFIIR.NE.N 
Ml  IR,1)=M< 
WIIR.l)  =CE 
Ml  IR,2)=L 
GO   TO    91 

92  Ml IR, 1 )=MC 
W{  IR,1)=DE 
K4=IXI IR)+ 
WC  =  XUR )-( 
Ml  IR,2)=K4 
W(  IR,2  )  =  W0 

91    CONTINUE 
C  FINAL    POIN 

DO  98  1=1, 
IF(I .NE.l) 
E2=X0/1.41 
E1=X(I J/1. 


HARACTERISTICS    OF    QUANTIZER    WITH    OFFSET 

F    STEPS 

IZED     STANDARD    DEVIATION    OF    X 

ATION    COEFICIENT    OF    X 

ATION    STEP 


EAL*8I A-H,0-Z) 

0134) ,X 134) , FUN t 10001 ),W I  18. 2), Ml  18 f 2) 

RAMETERS 

6)N 
2. 

PARAMETERS 
9 


0 
10 


-1.) )/VAR 

3 

+  1 

CHARACTERISTICS 
N 

+  1.)/2.0)*STP 
/2.0)*STP)+F/VAR 


INTEGRATICN 


=  1  ,M3 

LE.-5.0)     GO    TO    1000 

GE.+5.0)    GO    TO    1001 


1,N2 

1)    GO    TO    922 

5.0) /OEL+1.0 
K3-1 )*QEL+5.0 


2)    GO    TO    92 
IR-1,2)+1 

L-WO 


IR-1 ,2)+l 

L-WO 

5.0J/DEL+1 
K4-1 )*DEL+5.0 


T 
N 

GO    TO    99 
4 
414 
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GO  TO  100 
99  IF(I.NE.N)  GO  TO  101 
E2  =  E1 

E1=XF/1.414 
GO  TO  100 

101  E2  =  E1 

E1  =  XU)  /1.414 
100  R2  =  DERF(E1 )-DERF(E2  ) 
RR=0(I )**2*R2 
SS=SS+RR 
98  CCNTINUE 

RATT=SS/2.0 

RAT  =  DSORT{'RATT  ) 

WRITE(6.102)RATT,RAT 

102  FORMAT! 2F15.9) 
SY=-5.0 
WRITE(6,84)N1,N2 

84    FORM AT { •     • ,214) 

00    30    MM=1,21 

FORM    FUNCTION    PHI 

A=2.0*(  1.0-RH0**2) 

AA=DSORT(A) 

DO    2    J  =  1,L 

DO    1    I=1,N 

IF(I.NE.l)     GO    TO    5 

E2=( XO-RHO*SY) /AA 

E1=(X( I )-RHC*SY)/AA 

GO    TO    2  5 
5    IF(I.NE.N)    GO    TO    3 

E2  =  E1 

E1=(XF-RHC*SY)  /AA 

GO    TO    2  5 
3    E2  =  E1 

El=( X( I )-RHC*SY) /AA 
25    P,1  =  DERF(E1)-DERF(E2) 

R=Q(1)*R1 

S=S+R 

1  CONTINUE 

FUN( J)=(S/2.0)*DEXP((SY**2)/(-2.0) ) 

S=.0 

SY=SY+DEL 

2  CCNTINUE 

MAIN    INTEGRATION 

DO    7    I1=N1,N2 

IFCI1.NE.N1)    GO    TO    111 

M1=M(I1, 1) 

M2=M( 11,2) 

CORR=W(Ii,2)*FUN(M2) 

GO  TO  15 
111  IF<  I1.NE.N2 J  GO  TO  11 

M1=M( 11,1) 

M2=M(I1,2) 
■  CORR=W<  11,1  )*FUN(M1  ) 

GO  TO  15 
11  M1=M(I1,1) 

M2=M(I1,2) 

CORR=W( I  1,1)*FUN(M1)  +  W( I1,2)*FUN(M2) 
15  M3=M2-1 

DO  21  KC=M1,M3 

011=FUN(K0)+FUN(K0+1) 

E=E+011 
21  CCNTINUE 

RO=0(I1)*(DEL*B*0.5+CORR) 

FCQ=RQQ+RQ 
□  _  -\ 

7    CONTINUE 
ROOQ=RQQ*PP 
R000=(1. 0/RATT)*RQOG 
WRITE(6,40)RHO,RO0Q,R00O 
40     F0RMAT(3C25.9) 

WRITE!  7.204)  RCQCRCCQ 
204    FORMAT( 2F8.5) 
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RO0=0.0 

s=.o 

SY=-5.0 

IFCMM.NE.l )    GO    TO 

RHG=RHO+0.0499 

GO    TO    3  0 

67 

RFC=RH0+0.05 

30 

CONTINUE 

STOP 

END 

67 


C      DEGRADATION  FACTOR.   RC  FILTERING  AND  SAMPLING 
C      OFFSET  CASE 

IMPLICIT  RE£L*8(A-H,0-Z) 

REAL*4  F5,F6,F7,X3 

REAL  LABEL/AH     / 

REAL*8  ITITLEC12  J/'SWETT    Sll*'         ■/ 

REAL*8  TITLE11 JJ/10*'  •/ 

DIMENSION  Z(22)  ,F2C  22)  ,WH22),DELY1  22)  ,B(  22)  ,  SBC  22) 
1,FCC2,31),AP1C31),EC20),QC20),F5C31),F6C31)»F7C31), 
2X3C31)  ,F(2,31) 

GAUSS CE,VAR)=0. 3989423* ( DEXP (-C E/VAR)**2/2. ) )/VAR 

VAR=0.2 

FF=-0.5 

N=13 

AA=7. 2906748 

G=6. 24225363 

A=AA-G 

STP=I2/(N-1.0) ) 

DER=0.0 
C      CALCULATICN  OF  DERIVATIVE  OF  EXPECTED  VALUE 

DO  79  I=1,N 

0(I)  =  (  I-CN+1.0)/2.0)*STP 

EC  I)  =  C I-N/2.0)*STP+FF 

IFCI.NE.l)  GO  TO  2000 

E3=E( I) 

E1=0.0 

E2=GAUSSCE3,VAR) 

GO  TO  25U0 

2000  IFCI.NE.N)  GO  TO  2001 
E1=E2 

E2=0.0 

GO  TO  2500 

2001  E4=E(I) 
E1  =  E2 
E2=GAUSS (E4.VAR) 

2500  PPEL  =  0(  I  )*CE1-E2) 
79  CER=DER+PPEL 

.  WRITE(6.2531)A,DER 

2501  FORMATC  '  •  ,2F12.7) 
CO  98  1=1,21 

98  READ15.99) WI (I ) 

99  FORMAT CF 7.1) 
DO  1  1=1,21 
READ(5,2JF2<I ) 

1  CONTINUE 

2  F0RMATCF8.5) 
S=0.0 

CO    3    1=1.21 

Z(I)=S 

S=S+.05 

3  CONTINUE 

C      FIT  POLYNOMIAL 

CALL  LS0PL2{21,-11,Z,F2,WI ,Y,DELY, B, S3, TITLE) 

TRC=20.0 

BC1)=0.0 

XX=.01 

CO  35  1=1,31 

NN=TRC/XX-1.0 
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3332 
3333 


66 
35 


7 
987 


297 

888 

13 
31 

67 


RUN=0.0 
FRUN=0. 
DO  4  J  = 
T=XX*J 
X=DEXP( 
IF(X.LT 
PUN=B(  1 
1+B(6)*X 
2+B(ll)* 
GO  TO  3 
PUN  = (F2 
RUN=RUN 
RRUN=RR 
CONTINU 
A1=A^XX 
A2=XX*A 
A3=XX** 
F41=A1+ 
F42=A1+ 
FC(2,I ) 
F(2.  I)  = 
IFd.NE 
XX=XX+. 
GO  TO  3 

xx=xx+. 

CONT  INU 

XX  =  .l 

F(l,l)= 

YYY=1.0 

FC(1,1  ) 

X3( 1)=0 

DO  7  1  = 

YY=(XX/ 

Y4=(XX/ 

Y3=YY-( 

F(1,I)= 

FC(1,I  ) 

X3(I)=X 

XX=XX+. 

CONT  INU 

WRITE(6 

FORMAT( 

WRITEI6 

X=0.0 

DO  297 

RR=F(1, 

SS  =  (X*S 

APP=RR+ 

AP1(J)= 

X=X+.l 

CONTINU 

WRITE(6 

FORMAT! 

WRITE(6 

WRITEI6 

FORMAT( 

WRITEC6 

FORMAT! 

DO  67  K 

F5(K)  =  F 

F6<K)=F 

F7(K  i  =  A 

CALL  DR 

16,1, LAS 
CALL  DR 

16,1. LAS 
CALL  DR 

16,1, LAS 
STOP 
END 


0 

1,NN 


T) 

0.  1)  GO  TO  3332 

+B(2)*X+B(3  )*X**2+B(4)^X**3+B(  5)*X**4 

*5+e(7)*X**6+B(8 }*X**7+B(9)*X**8+B(10)*X**9 

**10+B( 12)*X**11 

33 

3)-F2( 1)  )*10.0*X 

PUN 

N+J*PUN 

2.0 

*RUN 

*AA*RRUN/TRC 

2- A  3 
2 

DSQRT(F41) /DER 
S0RT(F42)/DER 
1)     GO    TO    66 
9 


1.0 
-II. i 

=  CS0f 

.0 

2,31 

2.3) 

2.0) 

1  .0- 

DSOR' 

=  DS0f 

X 

1 

E 

,987 

•  1«  ) 

, 31)(X3t  I),F(1,  I), Ft  2, I) 1 1  =  1, 31) 


.O-DEXP(-TRC) J/TRC 
IRT(YYY) 


/DTANKMXX/2.0) 
/DSINHC XX/2»0) 
■DEXP(-TRC)  }*Y4**2/TRC 

:t(yy) 

)RT(  Y3) 


) 


J=l,31 
J  )**2 

TP**2)/24„0 
SS/VAR**2 

DSQRT(APP) 


,888) 

•  ',/////! 

,31)(X3(J),F(l,J),APKJ),J  =  lt31) 
,13) 

•  ',/////) 
,31)(X3U),FC(1,I),FC{2,I),I  =  1,31) 

•  »,3F12.7) 
=  1,30 
(1,K) 
<2,K) 
PKK) 

AW(30,X3,F6,l,0,LABEL,ITITLE,0.5tO.C,0,0,l,2,6, 
T) 

AW( 30.X3, F5,2, 0, LABEL, I TI TLE , 0. 5 , 0. C ,0  ,0 , 1 ,2 ,6, 
T) 

AW(33,X3,F7,3,4,LABEL,ITITLE,0.5,0.0,0,0,1,2,6, 
T  ) 
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DYNAMIC  RA 
IMPLICIT  R 
REAL*4  X3, 
REAL  LABEL 
REAL*8  ITI 
DIMENSION 
1F81250) , F9 
XO=-10.0E 
XF=10.0E  1 
DO  50  JJ=1 
READ(5,99) 

99    FORMAT (12) 
WRITE(6,83 

88  FORMAT( «  • 
S=.0 
VVAR=.l 
DO  8  L  =  l,2 
VAR=DSORT( 
DO  1  1=1, N 
X(I)=( I-N/ 

om=ii-(N 

IFU.NE.  1) 
E2=X0/1.41 
E1=X(I )/l. 
GO    TO    6 

5  IF(I.NE.N) 
E2  =  E1 

E1=XF/1.41 
GO    TO    6 

3    E2=E1 

E1=X(I  )/l. 

6  R1=DERF(£1 
R  =  Q(I )#*2# 
S  =  S+R 

1    CONTINUE 
QARR  =  WAR* 
Y(  JJ,L)=OA 
X3(LJ=VVAR 
PER(L)=( 1. 
VVAR=VVAR+ 
S=.0 
8    CONTINUE 
WRITE(6,82 
82    FORMATl'     « 
50    CONTINUE 
DO    46    K  =  l, 
F5(K)=Y(1, 
F6(K)  =  Y(2, 
F7(K)  =  Y(3, 
F8(K)=Y(4, 
F9(K)  =  Y(  5, 
46    F10(K)  =  Y(6 
CALL    DRAW( 

1,7,1, LAST) 
CALL    DRAW! 

1,7,1, LAST) 
CALL    DRAW I 

1,7,1,LAST) 
CALL    DRAWf 

1,7,1, LAST) 
CALL    DRA'w( 

I, 7,1, LAST) 
CALL    DRAW( 

16, 7,1, LAST 
STOP 
END 


NGE    OF    I 

EAL*8(  A- 

F5,F6,F7 

/4H  / 

TLE( 12)  / 

0(34), X( 

(250) ,F1 

10 

0 

,6 

N 

)N 
,12) 


50 
WAR) 

2.  )/VAR 
+1.J/2.J/VAR 

GO    TO    5 
4 
414 


F.    SAMPLING 
H,0-Z) 
,F8,F9,F10 

«  SWETT  •  ,11*'  •/ 

34),Y(6,250),X3(2  50),F5(25  0), 
0(25  0),PER(25  0),F6(250),F7(250) 


GO    TO    3 


414 

)-DERF(E2) 
Rl 


S 

RR 

0-QARR/(2.0*  WAR)  )*100. 
.1 


)(X3(L),Y(JJ,L),PER(L)t.L  =  l,250,3) 
•3F15.9) 


250 

K). 

K) 

K) 

K) 

K) 

,K) 

2  5  0,X3,F 


250,X3,F 

250,  X3,F 

250, X3,F 

250,X3,F 

250,X3,F 
) 


5,1, Of  LABEL, IT IT LE, 3. 0,5. 0,0, 0,2 ,2, 6 

6, 2, 0,  LAB  EL,  IT  ITLE, 3. 0,5. CO, 0,2, 2, 6 
7, 2,0,  LABEL,  ITI  TLE, 3.  0,5. CO, 0,2, 2, 6 
8, 2,0,  LABEL, I  TITLE  ,3. 0,5. CO, 0,2, 2, 6 
9, 2,0, LA BEL, ITITLE, 3. 0,5. 0,0, 0,2, 2, 6 
1C3,0,LABEL,ITITLE,3.C,5.C,0,0,2,2, 
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DYNAMIC  RANGE  OF  RC  FILTERING  AND  THEN  SAMPLING 
IMPLICIT  REAL*8( A-H.O-Z) 

RFAL*8  ITITLEC12  )/"  SWETT    'vll*'  '/ 

REAL  LABEL/4H     / 
REAL*4  X3, F5,F6,F7,F8,F9,F10 

DIMENSION  0(34) ,X(34)  ,Y ( 7 ,250 ) ,X3 (2501 , F5 (250 ) , 
1F8(250),F9(2  50),F10(250),PER(250),NN(8) ,VAN(5) 
2,F6(250) ,F7(250) 
SQ  =  5.0 
XF=10.0E  10 
XC=-10.0E  10 

READ(5,73MVAN(  I)  ,1=1,1) 
73  F0RMAT(F4.1) 

READ(5,71)  (NN<  I)  ,1=1,7) 

71  F0RMAT(7I2) 
1  =  1 
WRITE(6,76) VAN( I) 

76  FORMATl «  '  ,F4.1,//) 
VAR=VAN(  I) 
DO  70  JJ=1,6 
N=MN(JJ ) 
WRITE(6,72)N 

72  FORMAT(  •  '  ,12,//) 
ALPHA=1.0 

DC  19  K=l,110 

S=0.0 

CO  20  I  1=1, N 

X(I I)=( II-N/2.0) 

OUI)-(II-(N+1.0)/2.0) 

IFMI.NE.l  )  GO  TC  99 

E2=IX0+( VAR-ALPHA*VAR)*S0)/(1.414*VAR*ALPHA) 

E1  =  (X( I  I  )+{ VAR-ALPHA*VAR)*SQ)/( 1. 414* VAR* ALPHA) 

GO  TO  100 
99  IF(II.NE.N)  GO  TO  101 

E2  =  E1 

E1=(XF+(VAR-ALPHA*VAR)*SQ)/(1.414*VAR*ALPHA) 

GO  TO  100 
101  E2=E1 

E1=(X(II )+(VAR-ALPHA*VAR)*SQ)/ ( 1  .414*VAR* ALPHA) 
100  R1=DERF(E1)-DERF(E2) 

R2=0( II )*R1 

S=S+R2 
20    CONTINUE 

Y(JJ,K)=S/2.0 

X3(K)=ALPHA 

IFIK.EQ.1)     GO    TO    33 

PER(K)=I1.0-(Y ( JJ,K)/VAR)/((ALPHA-1.Q)*SQ) )*IC0.0 

GO    TO    34 

33  PER(l)=Q.O 

34  ALPHA=ALPHA+.l 
19    CONTINUE 

WRITE(6,75)(X3(L),Y(JJ,L),PER(L),L=1,110,3) 
75    FORMATC     •  ,3F15.7) 
70    CONTINUE 
9    CONTINUE 

DO   46    K=l,110 

F5(K)  =  Y(  1,K) 

F6(K)=Y(2,K) 

F7(K)=Y(3,K) 

F8(K  )  =  Y(  4,K) 

F9(K)=Y(5,K) 
46    F10(K)=YI6,K) 

CALL  DRAWC  110, X3,F 5, 1,0, LABEL, I  TITLE, 2. 0,2. 0,0, 0,2, 2,  6 
1  ,6,1, LAST) 

CALL  DR A U( 110, X3,F6, 2,0, LABEL, I TITLE, 2. 0,2. 5, 0,0, 2, 2, 6 
1,7,1, LAST) 

CALL  DRAW(110,X3,F7, 2,0, LABEL, IT ITLE,3.0, 5 .0, 0, 0, 2, 2, 6 
1,7,1, LAST) 

CALL  DRAW(110,X3,F8,2,0,LABEL,ITITLE,3.0,5.0,0,0,2,2,6 
1,7,1, LAST) 

CALL    DRAW( 110, X3,F 9, 2,0, LABEL, I  TITLE ,3. 0,5. 0,0, 0,2, 2, 6 
1,7,1, LAST) 
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CALL    DRAWC 1 1 0, X3  ,F 1 3 ,3 , 0 , LABEL , IT ITL E ,3 . 3 , 5 .0 , 0, 0 , 2, 2, 
16,7,1, LAST) 

STOP 
END 


DYNAMIC    RANG    OF    N.A.R.    RADIOMETER 

IMPLICIT    REAL*8(A-H,0-Z) 

REAL*4    X3,F5,F6,F7,F8,F9,F10 

REAL    LABEL/4H  / 

REAL*8     ITITLE(12)/« SWETT  «,11*«  •/ 

DIMENSION    Q(34) ,X(34) ,X3(250) ,F5(250) ,F6(250) ,F7(250) 

X0=-10.0E    1Q 

XF=10.JE    10 

READ(5,99)N 
99    F0RMAT(I2) 

WRITE! 6 • 88) N 
88    FORMAT! '     • ,12) 

S  =  .0 

VVAR=.l 

DO    8    L=l,250 

VAR=DSQRT(VVAR) 

DO    1    1=1, N 

X( I)=( I-N/2. )/VAR 

0(I)=( I-!N+1.) /2.J/VAR 

IF(I.NE.l)    GO    TO    5 

E2=X0/1.414 

E1=X(I J/1.414 

GO    TO    6 

5  IF( I.NE.N)    GO    TO    3 
E2=E1 
E1=XF/1.414 

GO    TO    6 
3    E2=E1 

E1=X(I) /1.414 

6  R1  =  DERF(E1 )-DERF(E2  ) 
R  =  QU  )**2*R1 

S  =  S+R 
1     CONTINUE 

OARR=VVAR*S 

F5(L)=0ARR 

X3(L)=VVAR 

VVAR=VVAR+.l 

S=.0 
8    CCNTINUE 

DO    87    K=l,180 

F6(K)=F5(K+40) 
87    F7(K)=F6(K)/F5 (K) 

WRITE(6,82){X3(L),F5(L),F6(L),F7(L),L=1,180,3) 
82     FORMAT!'     ' ,4F15.9) 

CALL    DRAW(183,X3,F5,1,3,LABEL,ITITLE,3.3,3.3,3,0,2,2,6 
1,7,,  1,  LAST) 

CALL    DRAW!  180.  X3,F6,2,0,  LABEL,  IT  ITLE, 3.0,  5. CO, 0,2, 2, 6 
1,7,1, LAST) 

CALL    DRAW(180,X3,F7,3,0,LABEL,ITITLE,3.0,5.0,0,0,2,2,6 
1  ,7,1 .LAST) 

STOP 

END 
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